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Abstract 

The splitting of invariant manifolds of whiskered (hyperbolic) tori with two frequencies in a nearly-integrable 
Hamiltonian system, whose hyperbolic part is given by a pendulum, is studied. We consider a torus with a fast 
frequency vector ui/y/e, with ui = (1, fi) where the frequency ratio fl is a quadratic irrational number. Applying the 
Poincare-Melnikov method, we carry out a careful study of the dominant harmonics of the Melnikov potential. This 
allows us to provide an asymptotic estimate for the maximal splitting distance, and show the existence of transverse 
homoclinic orbits to the whiskered tori with an asymptotic estimate for the transversality of the splitting. Both 
estimates are exponentially small in e, with the functions in the exponents being periodic with respect to lne, and 
can be explicitly constructed from the continued fraction of SI. In this way, we emphasize the strong dependence of 
our results on the arithmetic properties of SI. In particular, for quadratic ratios Si with a 1-periodic or 2-periodic 
continued fraction (called metallic and metallic-colored ratios respectively), we provide accurate upper and lower 
bounds for the splitting. The estimate for the maximal splitting distance is valid for all sufficiently small values 
of e, and the transversality can be established for a majority of values of e, excluding small intervals around some 
transition values where changes in the dominance of the harmonics take place, and bifurcations could occur. 

Keywords : splitting of separatrices, transverse homoclinic orbits, Melnikov integrals, quadratic frequency ratio. 


1 Introduction and setup 

1.1 Background and state of the art 

This paper is dedicated to the study of the exponentially small splitting of separatrices in a perturbed 3-degree- 
of-freedom Hamiltonian system, associated to a 2-dimensional whiskered torus (invariant hyperbolic torus) whose 
frequency ratio is an arbitrary quadratic irrational number (i.e. a real root of a quadratic polynomial with integer 
coefficients). 
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and the Russian Scientific Foundation grant 14-41-00044 at the Lobachevsky University of Nizhny Novgorod. The author MG has also 
been supported by the Swedish Knut and Alice Wallenberg Foundation grant 2013-0315. We also acknowledge the use of EIXAM, the UPC 
Applied Math cluster system for research computing (see http://www.mal.upc.edu/eixam/), and in particular Pau Roldan and Albert 
Granados for their support in the use of this cluster. 
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We start with an integrable Hamiltonian Hq having whiskered (hyperbolic) tori with coincident stable and unstable 
whiskers (invariant manifolds). We focus our attention on a torus, with a frequency vector of fast frequencies: 

= -^= , w = (l,fi), (1) 

v e 

whose frequency ratio H is a quadratic irrational number. If we consider a perturbed Hamiltonian H = Ho + pH i, 
where p is small, in general the whiskers do not coincide anymore. This phenomenon has got the name of splitting 
of separatrices, and is related to the non-integrability of the system and the existence of chaotic dynamics. If we 
assume, for the two involved parameters, a relation of the form p = e r for some r > 0, we have a problem of singular 
perturbation and in this case the splitting is exponentially small with respect to e. Our aim is to detect homoclinic 
orbits (i.e. intersections between the stable and unstable manifolds) associated to persistent whiskered tori, provide an 
asymptotic estimate for both the maximal splitting distance and its transversality , and show the dependence of such 
estimates on the arithmetic properties of the frequency ratio Q. 

To measure the splitting, it is very usual to apply the Poincare-Melnikov method , introduced by Poincare in [Poi90j 
and rediscovered much later by Melnikov and Arnold |Mel631 IArn64| . By considering a transverse section to the stable 
and unstable perturbed whiskers, one can consider a function A4(9), 6 £ T 2 , usually called splitting function , giving 
the vector distance, with values in R 2 , between the whiskers on this section, along the complementary directions. The 
method provides a first order approximation to this function, with respect to the parameter p, given by the (vector) 
Melnikov function M(9 ), defined by an integral (see for instance |Tre94 | ). We have 

M[9) = pM(9)+0(p 2 ), (2) 

and hence for p small enough the simple zeros of M(6) give rise to simple zeros of A i(9), i.e. to transverse intersections 
between the perturbed whiskers. In this way, we can obtain asymptotic estimates for both the maximal splitting 
distance as the maximum of the function \M(9)\, and for the transversality of the splitting, which can be measured 
by the minimal eigenvalue (in modulus) of the (2 x 2)-matrix DAd(0*), for any given zero 9*. 

An important fact, related to the Hamiltonian character of the system, is that both functions A 4(9) and M{9) are 
gradients of scalar functions |Eli941 IDGOOli : 

M(9) = V£(9), M(9) = VL(0). (3) 

Such scalar funtions are called splitting potential and Melnikov potential respectively. This means that there always 
exist homoclinic orbits, which correspond to critical points of C. 

As said before, the case of fast frequencies as in CD, with a perturbation of order p = e r , for a given r as small 
as possible, turns out to be a singular problem. The difficulty comes from the fact that the Melnikov function M(9) 
is exponentially small in e, and the Poincare-Melnikov method can be directly applied only if one assumes that p is 
exponentially small with respect to e. In order to validate the method in the case p = e r , one has to ensure that 
the error term is also exponentially small, and the Poincare-Melnikov approximation dominates it. To overcome such 
a difficulty in the study of the exponentially small splitting, it was introduced in |Laz03| the use of parameterizations 
of a complex strip of the whiskers (whose width is defined by the singularities of the unperturbed ones) by periodic 
analytic functions, together with flow-box coordinates. This tool was initially developed for the Chirikov standard 
map ILaz03l . and allowed several authors to validate the Poincare-Melnikov method for Hamiltonians with one and 
a half degrees of freedom (with 1 frequency) |HMS881 ISch891 IDS921IDS971 IGel97j and for area-preserving maps |DR98j . 

Later, those methods were extended to the case of whiskered tori with 2 frequencies. In this case, the arithmetic 
properties of the frequencies play an important role in the exponentially small asymptotic estimates of the splitting 
function, due to the presence of small divisors. This was first mentioned in |Loc90| . later detected in ISim94l . and 
rigorously proved in |DGJS97l for the quasi-periodically forced pendulum, assuming a polynomial perturbation in the 
coordinates associated to the pendulum. Recently, a more general (meromorphic) perturbation has been considered in 
IGS12j . It is worth mentioning that, in some cases, the Poincare-Melnikov method does not predict correctly the size 
of the splitting, as shown in |BFGS12l . where a Hamilton-Jacobi method is instead used. This method was previously 
used in SanOl . I'.MSO-l, KU’lloj . where exponentially small estimates for the transversality of the splitting were obtained, 
excluding some intervals of the perturbation parameter e. Similar results were obtained in lDG04l rDG03] . Moreover, 
the continuation of the transversality for all sufficiently small values of e was shown in IDG04I for the concrete case 
of the famous golden ratio Q = (y/b — l)/2, and in [DGG14c] for the case of the silver ratio = \[2 — 1, provided 


2 




































certain conditions on the phases of the perturbation are fulfilled. Otherwise, homoclinic bifurcations can occur, as 
studied, for instance, in ISVOlj for the Arnold’s example. Let us also mention that analogous estimates could also be 
obtained from a careful averaging out of the fast angular variables [PTOOl , at least concerning sharp upper bounds of 
the splitting. 

In general, in the quoted papers the frequency ratio is assumed to be a given concrete quadratic number (golden, 
silver). A generalization to some other concrete quadratic frequency ratios was considered in |DG03j . extending the 
asymptotic estimates for the maximal splitting distance, but without a satisfatory result concerning transversality. 
Recently, a parallel study of the cases of 2 and 3 frequencies has been considered in IDGG14ai (in the case of 
3 frequencies, with a frequency vector u> = (1,12,12 2 ), where Q is a concrete cubic irrational number), obtaining 
also exponentially small estimates for the maximal splitting distance. 

In this paper, we consider a 2-dimensional torus whose frequency ratio 12 in (O is a given quadratic irrational 
number. Our main objective is to develop a methodology, allowing us to obtain asymptotic estimates for both the 
maximal splitting distance and the transversality of the splitting, whose dependence on e is described by two piecewise- 
smooth functions denoted h\(e) and / 12 (e), respectively (see TheoremQ]). Such functions are both periodic with respect 
to lne, and their behavior depends strongly on the arithmetic properties of the frequency ratio 12. In particular, we 
show that the functions h\{e) and / 12 (e) can be constructed explicitly from the continued fraction of 12, and we can 
deduce some of their properties like the number of corners in each period (this can be seen as an indication of the 
complexity of the dependence on e of the splitting). Our goal is to show that our methods can be applied to an 
arbitrary quadratic ratio, and hence the results on the splitting distance and transversality generalize the ones of 
)PG03| DG04 . DGG14a ; DGG14cI . Although we do not study here the continuation of the transversality for all 
values of £ —X 0, we stress that this could be carried out by means of a specific study in each case, as done in 
jSVOli .1K10 1. lDGG14c i for some concrete (golden, silver) ratios. 

We point out that the periodicity in In e of the functions h\ (e) and /12 (&) comes directly from the special properties 
of the continued fractions of quadratic numbers and cannot be satisfied in other cases (see |DGG14b] . where the case 
of numbers of constant type is considered). 

As we describe in Section [ITSl the dependence on e of the function /i 2 (£), associated to the asymptotic estimates 
for the transversality of the splitting, is more cumbersome than for the function h\(e), associated to the maximal 
splitting distance. We stress here that, for some purposes, it is not necessary to establish the transversality of the 
splitting, and it can be enough to provide estimates of the maximal splitting distance. Indeed, such estimates imply 
the existence of splitting between the invariant manifolds, which provides a strong indication of the non-integrability 
of the system near the given torus, and opens the door to the application of topological methods [GR03 , GL061 for 
the study of Arnold diffusion in such systems. 


1.2 Setup 

Here we describe the nearly-integrable Hamiltonian system under consideration. In particular, we study a singular or 
weakly hyperbolic (a priori stable) Hamiltonian with 3 degrees of freedom possessing a 2-dimensional whiskered torus 
with fast frequencies. In canonical coordinates ( x , y,ip,I) €E TxlxT 2 xI 2 , with the symplectic form da;Ady + d</?Ad/, 
the Hamiltonian is defined by 


H(x,y,ip,I) = H 0 (x,y,I) + yH^x, <p), (4) 

Ho(x,y,I) = (w e ,7) + ^(A 1,1) + y +cosx - 1, (5) 

Hi{x,<p) = h(x)f(<p). (6) 

Our system has two parameters e > 0 and y, linked by a relation y = s r , r > 0 (the smaller r the better). Thus, if we 
consider e as the unique parameter, we have a singular problem for e —X 0. See |DG01| for a discussion about singular 
and regular problems. 

Recall that we are assuming a vector of fast frequencies uj e = oj/y/e as given in dTJ) , with a frequency vector 
u> = (1,12) whose frequency ratio 12 is a quadratic irrational number ; we assume without loss of generality that 
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0 < fl < 1. It is a well-known property (and we prove it in Section HOI see also ILan95l §11.2]) that any vector with 
quadratic ratio satisfies a Diophantine condition 

\(k,u)\>^ 7 r Viez 2 \{ 0 }, (7) 

with some 7 > 0. We also assume in ([5]) that A is a symmetric (2 x 2 )-matrix, such that Hq satisfies the condition of 
isoenergetic nondegeneracy 

det ( ^ o)*°- (8) 

For the perturbation Hi in ([ 6 ]). we deal with the following analytic periodic functions, 

h(x) = cosx, /(</?) = e~ p \ k \ cos((fc, (p) — (Jk), with 07 £ T, (9) 

k£Z 

where we introduce, in order to avoid repetitions in the Fourier series, the set 

Z = {k = (ki,k 2 ) £ Z 2 : k 2 > 0 or ( k 2 = 0, k\ > 0)}. (10) 

The constant p > 0 gives the complex width of analyticity of the function f(p). Concerning the phases 07 , they can 
be chosen arbitrarily for the purpose of this paper. 

To justify the form of the perturbation Hi chosen in (ED and we stress that it makes easier the explicit 
computation of the Melnikov potential, which is necessary in order to show that it dominates the error term in ([2]) , 
and therefore to establish the existence of splitting. Moreover, the fact that all coefficients fk = e~ p ^, in the Fourier 
expansion with respect to <p, are nonzero and have an exponential decay, ensures that the study of the dominant 
harmonics of the Melnikov potential can be carried out directly from the arithmetic properties of the frequency 
vector w (see Section (Sj. Since our method is completely constructive, a perturbation with another kind of concrete 
harmonics fk could also be considered (like fk = \k\ m e~ p ^), simply at the cost of more cumbersome computations 
in order to determine the dominant harmonics of the Melnikov potential. 

It is worth reminding that the Hamiltonian defined in (J4][9D is paradigmatic, since it is a generalization of the 
famous Arnold’s example (introduced in |Arn64| to illustrate the transition chain mechanism in Arnold diffusion). It 
provides a model for the behavior of a nearly-integrable Hamiltonian system near a single resonance (see IDG01I for 
a motivation), and has often been considered in the literature (see for instance [GGM99al ILMS031IDGS04] '). We also 
mention that a perturbation with an exponential decay as the function f(tp) in d9| has also been considered (see for 
instance IPTOOj l. In the present paper, our aim is to emphasize the dependence of the splitting, and its transversality, 
on the arithmetic properties of the frequency vector w. 

Let us describe the invariant tori and whiskers, as well as the splitting and Melnikov functions. First, it is clear that 
the unperturbed system H 0 (that corresponds to p, = 0) consists of the pendulum given by P(x,y) = y 2 /2 + cosa;— 1, 
and 2 rotors with fast frequencies: ip = to e + A I, 1 = 0. The pendulum has a hyperbolic equilibrium at the origin, 
with separatrices that correspond to the curves given by P(x, y) = 0. We parameterize the upper separatrix of the 
pendulum as (xo(s), yo{s)), s € R, where 

2 

Xq(s ) = 4arctane s , Vo(s) = -:— • 

coshs 

Then, the lower separatrix has the parametrization (xo(—s), — yo{— s)). For the rotors system (</?,/), 
are I = Iq, p = po + ^(<W + A/o). Consequently, the Hamiltonian Hq has a 2-parameter family of 
whiskered tori: in coordinates (x,y, (/?, I), each torus can be parameterized as 

Ti 0 : (0,0,0, J 0 ), 9 £ T 2 , 

and the inner dynamics on each torus is 9 = oj e + A/o. Each invariant torus has a homoclinic whisker , 

3-dimensional stable and unstable invariant manifolds, which can be parameterized as 

W / 0 : (x 0 (s),yo(s),9,I 0 ), s € R, 9 £ T 2 , ( 11 ) 


the solutions 
2 -dimensional 


i.e. coincident 
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with the inner dynamics given by s = 1, 6 = oj e + AIq. 

In fact, the collection of the whiskered tori for all values of Iq is a 4-dimensional normally hyperbolic invariant 
manifold , parameterized by (0,1) e T 2 x M 2 . This manifold has a 5-dimensional homoclinic manifold, which can be 
parameterized by ( s,0,I ), with inner dynamics s = 1, 0 = oj £ + A I, 1 = 0. We stress that this approach is usually 
considered in the study of Arnold diffusion (see for instance [DLS061 1. 

Among the family of whiskered tori and homoclinic whiskers, we will focus our attention on the torus 7o, whose 
frequency vector is uj e as in m, and its associated homoclinic whisker Wo. 

When adding the perturbation fiH\, for /i ^ 0 small enough the hyperbolic KAM theorem can be applied (see for 
instance |NieOO] l thanks to the Diophantine condition (0 and to the isoenergetic nondegeneracy ([8]). For /i small 
enough, the whiskered torus persists with some shift and deformation, as a perturbed torus T = as well as its 

local whiskers Wi oc = W^ (a precise statement can be found, for instance, in [DGS041 Th. 1]). 

The local whiskers can be extended along the flow, but in general for p ^ 0 the (global) whiskers do not coincide 
anymore, and one expects the existence of splitting between the (3-dimensional) stable and unstable whiskers, denoted 
W s = y\i s dv) and W u = W u ’^ respectively. Using flow-box coordinates (see |DGS04j . where the n-dimensional case is 
considered) in a neighbourhood containing a piece of both whiskers (away from the invariant torus), one can introduce 
parameterizations of the perturbed whiskers, with parameters (s, 0) inherited from the unperturbed whisker (1111) . and 
the inner dynamics 

s = 1, 0 = uj e . 

Then, the distance between the stable whisker W s and the unstable whisker W u can be measured by comparing such 
parameterizations along the complementary directions. The number of such directions is 3 but, due to the energy 
conservation, it is enough to consider 2 directions, say the ones related to the action coordinates. In this way, one can 
introduce a (vector) splitting function , with values in M 2 , as the difference of the parameterizations ff s ’ u (s,0) of (the 
action components of) the perturbed whiskers W s and W u . Initially this splitting function depends on ( s,6 ), but it 
can be restricted to a transverse section by considering a fixed s, say s = 0, and we can define as in [DG00| §5.2] the 
splitting function 

M(0) := J u (0,0) — J s (0,0), 6> € T 2 . (12) 

As said in ©, this function turns out to be the gradient of the (scalar) splitting potential C(6) (see IDGOOI IEli94l 1. 
Notice that the nondegenerate critical points of C correspond to simple zeros of A4 and give rise to transverse homoclinic 
orbits to the whiskered torus. 


Applying the Poincare-Melnikov method, the first order approximation © of the splitting function is given by 
the (vector) Melnikov function M(0), which is the gradient of the Melnikov potential: M(6) = X/L(6). The latter 
one can be defined as an integral: we consider any homoclinic trajectory of the unperturbed homoclinic whisker Wo 
in HUb starting on the section s = 0, and the trajectory on the torus 7o to which it is asymptotic as t — > ±oo, and 
we substract the values of the perturbation H\ on the two trajectories. This gives an absolutely convergent integral, 
depending on the initial phase 0 £ T 2 of the considered trajectories: 


m 


/ oo 

[Hi(x 0 (t),6 + tu e ) - H 1 (0,0 + tco e )\dt 

-oo 

/ oo 

[h(x 0 (t)) - h(0)]f(0 + tu e ) d t, 

-oo 


(13) 


where we have taken into account the specific form © of the perturbation. 


Our choice of the pendulum P(x, y ) = y 2 /2 + cosx — 1 in ©, whose separatrix has simple poles, makes it possible 
to use the method of residues in order to compute the coefficients L& of the Fourier expansion of the Melnikov 
potential L(8). Such coefficients turn out to be exponentially small in e (see their expression in Section [30). For each 
value of e only a finite number of dominant harmonics are relevant to find the nondegenerate critical points of L(0), 
i.e. the simple zeros of the Melnikov function M(0). Due to the exponential decay of the Fourier coefficients of f(ip) 
in ©, it is not hard to study such dominance and its dependence on e. 


In order to give asymptotic estimates for both the maximal splitting distance and the transversality of the ho¬ 
moclinic orbits, the estimates obtained for the Melnikov function M(6) have to be validated also for the splitting 
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function A4(8). The difficulty in the application of the Poincare-Melnikov approximation ([2]). due to the exponential 
smallness in e of the function M (9) in our singular case p = e r , can be solved by obtaining upper bounds (on a complex 
domain) for the error term in ([2]). showing that, if r > r* with a suitable r*, its Fourier coefficients are dominated by 
the coefficients of M(9) (see also [DGS04) ). 

We stress that our approach can also be directly applied to other classical 1-degree-of-freedom Hamiltonians 
P{x,y) = y 2 /2 + V(x), with a potential V(x) having a unique nondegenerate maximum, although the use of residues 
becomes more cumbersome when the complex parameterization of the separatrix has poles of higher orders (see some 
examples in [DS97] ). 


1.3 Main result 

For the Hamiltonian system mm with the 2 parameters linked by p = e r , r > r* (with some suitable r*), and 
a frequency vector u> = (1,12) with a quadratic ratio 12, our main result provides exponentially small asymptotic 
estimates for some measures of the splitting. On one hand, we obtain an asymptotic estimate for the maximal 
distance of splitting, given in terms of the maximum size in modulus of the splitting function A4(9) = V£(0), and 
this estimate is valid for all e sufficiently small. On the other hand, we show the existence of transverse homoclinic 
orbits, given as simple zeros 8* of A4(6) (or, equivalently, as nondegenerate critical points of C{6)), and we obtain an 
asymptotic estimate for the transversality of the homoclinic orbits, measured by the minimal eigenvalue (in modulus) 
of the matrix D A4(8*) = D 2 C(0*), at each zero of A4(0). This result on transversality is valid for “almost all” e 
sufficiently small, since we have to exclude a small neighborhood of a finite number of geometric sequences where 
homoclinic bifurcations could take place. 

With our approach, the Poincare-Melnikov method can be validated for an exponent r > r* with r* = 3, although 
a lower value of r* can be given in some particular cases (see remark [T] after Theorem [T]). However, such values of r* 
are not optimal and could be improved using other methods, like the parametrization of the whiskers as solutions of 
Hamilton-Jacobi equation (see for instance [LMS03 |. iBFGS12] l. In this paper, the emphasis is put in the generalization 
of the estimates to the case of an arbitrary quadratic frequency ratio fl, rather than in the improvement of the value 
of r*. 

Due to the form of f(<p) in ©, the Melnikov potential L{8) is readily presented in its Fourier series (see Section l3Jll . 
with coefficients Lk = Lk(e) which are exponentially small in e. We use this expansion of L{8) in order to detect its 
dominant harmonics for every given e. A careful control of the error term in d2|) ensures that the dominant harmonics 
of L(8) correspond to the dominant harmonics of the splitting potential C{9). Such a dominance is also valid for the 
splitting function M(9), since the size of their Fourier coefficients Mk (vector) and Ck (scalar) is directly related: 
\Mk\ = |fc| Ck, k € Z (recall the definition (ITOT) 1. 

As shown in Section [H in order to obtain an asymptotic estimate for the maximal distance of splitting, it is enough 
to consider the first dominant harmonic, given by some vector in Z which depends on the perturbation parameter: 
k = S i(e). Using estimates for this dominant harmonic Cs 1 as well as for all the remaining harmonics, we show that 
the dominant harmonic is large enough to ensure that it provides an approximation to the maximum size of the whole 
splitting function (see also [DGG14at IDGG14b| b On the other hand, to show the transversality of the splitting, it 
is necessary to consider at least two dominant harmonics in order to prove the nondegeneracy of the critical points 
of the splitting potential (see also [DG031100041 1. For most values of the parameter e, it is enough to consider the 
two “essential” dominant harmonics Cs 1 and Cg 2 (he. the two most dominant harmonics whose associated vectors 
Si (e), S 2 (e) € Z are linearly independent, see Section f?~?l) . and the transversality is then directly established. 

However, one has to consider at least three harmonics for e near to some “transition values” e, introduced below as 
the values at which a change in the second essential dominant harmonic occurs and, consequently, the second and some 
subsequent harmonics have similar sizes. Such transition values turn out to be corners of the function h 2 {s), related 
to the size of the second dominant harmonic (see the theorem below), and are given by a finite number of geometric 
sequences. The study of the transversality for e close to a transition value, which is not considered in this paper, 
requires to carry out a specific study that depends strongly on the quadratic frequency ratio 12, and on the concrete 
perturbation considered in ([9|l . In some cases, the continuation of the transversality for all sufficiently small values 
e —► 0 can be established under a certain condition on the phases er^ in ([9]>. as done in |DG04| and |DGG14c1 for the 
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golden and silver ratios, respectively. Otherwise, homoclinic bifurcations can occur when e goes accross a transition 
value (see for instance jSVOl] ). 


The dependence on e of the size of the splitting and its transversality is closely related to the arithmetic properties 
of the frequency vector ui = (l,f2), since the integer vectors k £ Z associated to the dominant harmonics can be 
found, for any e, among the main quasi-resonances of the vector w, i.e. the vectors k giving the “least” small divisors 
|(fc, w)| (relatively to the size of |fc|). In Section O we develop a methodology for a complete study of the resonant 
properties of vectors with a quadratic ratio, which is one of the main goals of this paper. This methodology relies 
in the classification, established in |DG03| for any vector with a quadratic ratio f2, of the integer vectors k into 
“resonant sequences” (see also Sections 12.21 and [2.31 for definitions). Among them, the sequence of primary resonances 
corresponds to the vectors k which fit best the Diophantine condition m, and the vectors k belonging to the remaining 
sequences are called secondary resonances. 

As particular cases, for the golden ratio = (\/5 — l)/2 the primary resonances can be described in terms of the 
Fibonacci numbers: k = (—F n -i,F n ) (see for instance lDG04t ). and in the case of the silver ratio Q = y/2 — 1 the 
primary resonances are given in terms of the Pell numbers (see |DGG14c] j. which play an analogous role. In general, 
for a given quadratic ratio f 1 the sequence of primary resonances, as well as the remaining resonant sequences, can 
be determined from the continued fraction of f2, which is eventually periodic, i.e. periodic starting at some element 
(see Section [2.ID . We can construct, from the continued fraction, a unimodular matrix T (i.e. with integer entries and 
determinant ±1), having ui as an eigenvector with the associated eigenvalue 

A = A(ft) > 1 

(see Proposition U for an explicit construction). Then, the iteration of the matrix (T -1 ) t from an initial ( “primitive”) 
vector allows us to generate any resonant sequence (see the definition (1211) I. 

Next, we establish the main result of this work, providing two types of exponentially small asymptotic estimates 
for the splitting, as e —> 0, related to the maximal distance of splitting and its transversality. The first one is given by 
the maximum of \A4(9)\, 9 £ T 2 . On the other hand, we show that for most values of e the function A 4(9) has simple 
zeros, which correspond to transverse homoclinic orbits, and for each zero 9* we give an estimate for the minimum 
eigenvalue (in modulus) of the matrix DA4(9*), as a measure for the transversality of the splitting. This generalizes 
the results of [DG031IDG04) . 

We stress that the dependence on e of both asymptotic estimates is given by the exponent 1/4, and by the functions 
h\(e) and /12(c), which are periodic with respect to lne and piecewise-smooth and, consequently, have a finite number 
of corners (i.e. jump discontinuities of the derivative) in each period. Some examples are shown in Figures |T] [2] (where 
a logarithmic scale for e is used). The oscillatory behavior of the functions h\{e) and /12(e) depends strongly on the 
arithmetic properties of f l and, in fact, both functions can be explicitly constructed from the continued fraction of Q 
(see Section 13.2D . Below, in two additional results we establish more accurately the behavior of the functions hi(e) 
and /12(e) in the simplest cases of 1-periodic and 2-periodic continued fractions. 

For positive quantities, we use the notation / ~ g if we can bound c\g < / < C 2 g with constants Ci,C 2 > 0 
not depending on e, /i. 


Theorem 1 ( main result) Assume the conditions described for the Hamiltonian with a quadratic frequency 

ratio fl, that e is small enough and that p = e r , r > 3. Then, for the splitting function A i(9) we have: 


(a) max \M(9)\ ~ —r-^exp 
ep t 2 n e 1 / 2 


Co /11 (e)) 
e 1 / 4 j’ 


(b) the number of zeros 9* of A4 (9) is 4k with u{e) > 1 integer, and they are all simple, for any e except for a small 
neighborhood of some transition values e, belonging to a finite number of geometric sequences; 


(c) at each zero 9* of M.{9), the minimal eigenvalue of T)A4(9*) satisfies 


| m 


pe 1 ^ 4 exp 


C 0 h 2 (e) \ 

e 1 / 4 J 
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The functions h\{e) and /12(e), defined in m are piecewise-smooth and 4 In A -periodic in lne, with A = A (f 2 ) as 
given in Proposition [7} In each period, the function /7.1(e) has at least 1 corner and /12(e) has at least 2 corners. They 
satisfy for e > 0 the following bounds: 

min/i 1 (e) = 1, max/i^e) < J 1 , max/i 2 (e) < J 2 , /ii(e) < /12(e), 

with the constants 

Ji = J 2 = J 2 (n):=^A+i). (14) 

The corners of hi(e) are exactly the points S such that /71(e) = /12(e). The corners of /12(e) see the same points i, 
and the points e where the results of (b-c) do not apply. The (integer) function k(e) is piecewise-constant and 41 nA- 
periodic in lne, eventually with discontinuities at the transition points e. On the other hand, Co = Cq( 0, p) is a 
positive constant defined in m • 




Figure 1: Dependence on e of the functions in the exponents, for the metallic ratio Q. = [3] (the bronze ratio), using a loga¬ 
rithmic scale for e, 

(a) graphs of the functions g* s ( qn ){e), associated to essential (the solid ones) and non-essential (the dashed ones) reso¬ 
nances s(q,n), the red ones correspond to the primary functions g n (s) (see Section \3.1\ ) : 

(b) graphs of the functions hi{e) and /12(e)- 




h^e) 


Figure 2: Graphs of the functions hi{e) and /i 2 (e) for two metallic-colored ratios, 


(a) = [1,3] (a golden-colored ratio); (b) fZ = [2,3] (a silver-colored ratio). 





















































Remarks. 


1. If the function h{x) in ([21) is replaced by h(x) = cos a; — 1, then the results of Theorem Q] are valid for p = e r 
with r > 2 (instead of r > 3). The details of this improvement are not given here, since they work exactly as 

in rurm 

2. As a consequence of the theorem, replacing h\(e) by its upper bound J\ > 0, we get the following lower bound 
for the maximal splitting distance: 

where c is a constant. This may be enough, if our aim is only to prove the existence of splitting of separatrices, 
without giving an accurate description for it. 

3. The results of Theorem [Q can be partially generalized if the frequency ratio fl is a non-quadratic number of 
constant type , i.e. whose continued fraction has bounded entries, but it is not periodic. The numbers of constant 
type are exactly the ones such that w = (1, fl) satisfies a Diophantine condition with the minimal exponent, as 
in 0. This case has been considered in [DGG14bl . where a function analogous to 6i(e), providing an asymptotic 
estimate for the maximal splitting distance, is defined. In this case, this function is bounded but it is no longer 
periodic in lne. 


For the simplest cases of continued fractions, we can obtain more accurate information on the functions hi{e) 
and 62(e). As we show in Section 12.11 we can restrict ourselves to the case of a purely periodic continued fraction, 
= [ai,..., a m ] (we assume that 0 < O < 1, see Section l2Jl for the notation). In particular, we consider the following 
two cases: 


• the metallic ratios: = [a] 


s/a 2 + 4 - a 
2 


with a > 1; 


.—ri V a2 b 2 + 4 ab — ab 

• the metallic-colored ratios: il = a,o =-, with 1 < a < b 

2 a 

(for the metallic-colored ratios, it is not necessary to consider the case a > b, since [a, b] = [a, b, a]). 


The metallic ratios, which are limits of the sequence of quocients of consecutive terms of generalized Fibonacci 
sequences, have often been considered (see for instance |Spi99[ IFP07j l. As some particular cases, we mention the 
golden , silver and bronze ratios: [1], [2] and [3], respectively. The metallic-colored ratios can be subdivided in 
several classes, such as: 

• the golden-colored ratios: = [ 1 , b ], with b > 2; 

• the silver-colored ratios: Ul = [2,6], with b > 3; 

• the bronze-colored ratios , etc. 

For such types of ratios, we next provide descriptions of the functions h\{e) and 62(e). Additionally, in part (b) 
we include a statement concerning the exact number of critical points of the splitting function Ai ( 9 ) , for the case of 
metallic ratios. Such results come from an accurate analysis of the first and second essential dominant harmonics of 
the Melnikov potential, studying whether they are both given by primary resonances for any e, or they can be given 
by secondary resonances for some intervals of £. We point out that a rigorous analysis of the role of the secondary 
resonances becomes too cumbersome in some cases. For this reason, although we give rigorous proofs for some results 
(the statements of Theorem 0, we provide numerical evidence for other ones (given below as “Numerical Result [3]’ ) 
after checking them with intensive computations carried out for a large number of frequency ratios (see the proofs and 
justifications in Section EPI . 


Theorem 2 Under the conditions of Theorem[ we have: 

(a) If the frequency ratio f l is metallic, the function h\ (e) has exactly 1 corner £ in each period, satisfies max hi (e) = 
hi (e) = Ji, and the distance between consecutive corners is exactly 4 In A. 


9 
















(b) If the frequency ratio fl is golden-colored, the function /12(e) has at least 3 corners in each period, and satisfies 
min/12(e) = J\, max/12(e) < J2- 

(c) If the frequency ratio is metallic-colored but not golden-colored, the function /11(e) has at least 2 corners in 
each period, and satisfies max/11(e) < J\. 


Numerical Result 3 Under the conditions of Theorem ]! J we have: 

(a) If the frequency ratio f l is golden-colored, the function h\{e) has exactly 1 corner e in each period, satisfies 

max/11(e) = /11(e) = J\, and the distance between consecutive corners is exactly 4 In A. 0 

(b) If the frequency ratio 0 is metallic, the function /12(e) has exactly 2 corners e, e in each period, satisfies 

min /12(e) = /12(e) = J\, max/12(e) = /12(e) = J2, and the distance between consecutive corners is exactly 2 In A. 

Moreover.the number of zeros 0 * of A 4 ( 9 ) is exactly f, for any e except for a small neighborhood of the transition 
values e. n 


As said before, the functions /11(e) and /12(e) can be defined explicitly for any quadratic ratio f2, from its continued 
fraction (see Section [3]). Such functions have piecewise expressions, which are simple in the case of a metallic ratio, 
but in general they can be very complicated, depending on the number of their corners in each period. Nevertheless, 
we stress that numerical justifications are required for results concerning infinite families of ratios, such as the metallic 
or the golden-colored ones. Instead, for particular frequency ratios the results can be rigorously established. 

Organization of the paper. We start in Section [2] with studying the arithmetic properties of frequency vectors 
w = (1,12) with a quadratic ratio f l. Such properties are closely related to the continued fraction of 0 /Section [2.11) . 
which allows us to construct the iteration matrices allowing us to study the resonant properties of the vector lo 
/Sections 12.21 and liPl) . and to provide accurate results for the cases of metallic and metallic-colored ratios /Section l2l4ll . 
mainly considered in this paper. Next, in Section [3] we find an asymptotic estimate for the first and second dominant 
harmonics of the splitting potential, which allows us to define the functions hi(e) and /12(e) and study their general 
properties /Sections 13.11 and 13.21) . as well as the specific properties for 1-periodic and 2-periodic continued fractions 
/Section 13.31) . considered in Theorem [2] and Numerical Result [3] Finally, in Section 2] we provide rigorous bounds of 
the remaining harmonics allowing us to obtain asymptotic estimates for both the maximal splitting distance and the 
transversality of the splitting, as established in Theorem [T| 

Finally, we introduce some notations that we use in this paper. For positive quantities, we write / A g if we 
can bound / < eg with some constant c not depending on e and g. In this way, we can write / ~<?if<7^/^g. 
On the other hand, when comparing positive sequences a n , b n we use an expression like “a n « b n as n —> 00 ” if 
lim ( a n /b n ) = 1, and also “ a n < b n as n —> 00” if limsup(a n /6„) < 1. 

n—y 00 r>.—^00 


2 Vectors with quadratic ratio 

2.1 Continued fractions of quadratic numbers 

It is well-known that any irrational number 0 < 12 < 1 has an infinite continued fraction 

0= [ai,a 2 ,a 3 ,...] =-*—j-. aj e Z + , j > 1 

ai H-j- 

<12 H-;- 

a 3 H- 

(notice that the integer part is a 0 = 0, hence we have removed the entry ‘0;’ from the notation). Its entries aj are 

T) ■ 

called the partial quotients of the continued fraction. It is also well-known that the rational numbers — = [ai,..., a^], 

lr The result of part (a) has been checked numerically for golden-colored ratios Q = [l,b], 2 < b < 10 6 . 

2 The results of part (b) have been checked numerically for metallic ratios = [a], 1 < a < 10 4 . 
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j > 1, called the (principal) convergents of f l, provide successive best rational approximations to Q. Thus, if we 
consider the “vector convergents” w(j) := ( qj,Pj ), we obtain approximations to the direction of the vector w = ( 1 , U) 
(see, for instance, |Sch80l and Lan95l as general references on continued fractions). 

The convergents of a continued fraction are usually computed from the standard recurrences 

q -1 = 0 , q 0 = 1 , qj = djqj -1 + ^_ 2 , 

(15) 

P -i = l, Po = a o = 0, pj = ajPj-i + pj_2, j> 1- 

Alternatively, we can compute them in terms of products of unimodular matrices |DGG14bl Prop. 1], 

=A 1 ---A j , where Ai = T(ai) ^ ^ (16) 

If we consider the first column, we can write w(j) = A\ • • • Aj'ui(O). 

An important tool in the study of continued fractions is the Gauss map g : (0,1) —> [0,1), defined as g(x) = 

where {•} stands for the fractional part of any real number. This map acts on a given continued fraction by removing 
the first entry: for LI = [ai, a 2 , 03 ,...], we have g(Ll) = [a 2 , 03 ,...]. We consider, for a given number LI £ (0,1), the 
sequence (xj) defined by 

x 0 = Ll Xj = g(xj- 1 ), j > 1, (17) 

which satisfies that Xj 0 for any j if LI is irrational. It is clear that Xj = [ aj+ 1 , aj- j_ 2 ,... ] for any j. 

In our case of a quadratic irrational number fl, it is well-known that the continued fraction is eventually periodic , 
i.e. periodic starting at some partial quotient. For an m-periodic continued fraction, we use the notation 

II — [bl, . . . ,b r , Ctl, . . . , dm ] ■ 

In fact, as we see below we can restrict ourselves to the numbers with purely periodic continued fractions, i.e. periodic 
starting at the first partial quotient: LI = [ aq,..., a m ]. It is easy to relate such properties with the sequence (xj) 
defined by the Gauss map: the continued fraction of LI is eventually periodic (hence, Ll is quadratic) if and only if 

x r +m = x r for some r > 0 , m > 1 , and it is purely periodic if and only if x m = xo for some m > 1 . 

In the following proposition, which plays an essential role in the results of this paper, we see that for any given 
vector a; = (1, Ll) with a quadratic ratio f2, there exists a unimodular matrix T = T(f2) having ui as an eigenvector with 
the associated eigenvalue A = A(fl) > 1. We show how we can construct both T and A, directly from the continued 

fraction of Q. Additionally, we show that applying the matrix T to a convergent w(j) we get the convergent w(j + m). 



( qj qj -1 
V Pj Pj -1 


Proposition 4 

(a) Let Li £ (0,1) be a quadratic irrational number with a purely periodic continued fraction: = [ai,..., a m ], 

and consider the matrices Aj = T~(a,j) as in 1161 ). Then, the matrix T = A\ ■ ■ ■ A m is unimodular, and has 

u: = (l,f2) as eigenvector with eigenvalue A = - > 1, where ( Xj ) is the sequence defined by 1171). 

X0Xl---Xm~l 

Moreover, for the convergents w(j) of f2 we have that Tw{j) = w(j + m) for any j > 0. 

(b) Let fl be a quadratic irrational number with a non-purely periodic continued fraction: fl = [ 61 ,..., b r , fl], with Q 
as in (a), and consider the matrices Bj = T(bj), and S = Bi-- B r . Then, the matrix T = STS~ l is 
unimodular, and has a) = (l,f2) as eigenvector with eigenvalue A as in (a). Moreover, for the convergents w{j) 
of fl we have that Tw(j) = w(j + m) for any j > r. 


Proof. Using the construction of the sequence (xj) associated to 0, we see that 


deduce the equality 


( 1 ' 

< 

\ 

A 

II 

' 1 ) 

, n > 1 

V Xj -1 J 


v X J J 



- = a 7 - + Xj, and we easily 

Xj -1 
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Iterating this equality for j = 1 ,,m and using that xq = fl = x m , we obtain 

= x 0 xi ■ ■ ■ x m -i AiA 2 ■■■A, 

which proves that Tw = Aw, and it is clear that T is unimodular. To complete part (a), using (fl6l) and the periodicity 
of the continued fraction we have 

Tw(j) = A x ■ ■ ■ A m A i • • • Ajw( 0) = Ai ■ • • A j+m w( 0) = w(j + m). 




With similar arguments we prove part (b). Indeed, using the sequence (xj) associated to fl, we see that 

^ = xoxi ■ ■ ■ x r -i B\B 2 ■ ■ ■ B r 

which says that the matrix S provides a unimodular linear change between the directions of the vectors w and w. We 
deduce that Tw = Aw. On the other hand, the matrix S also provides a relation between their respective convergents. 
Indeed, using (fTBl) we see that, for j > r, 

w(j) = Bi ■ ■ ■ B r A\ ■ ■ ■ Aj_ r w( 0) = Sw(j — r ) 

(notice that w(0) = w;(0) = (1,0)). Then, using (a) we deduce that 

Tw(j) = STw(j — r) = Sw(j) = w(j + r). 




□ 


Remarks. 

1. In what concerns the contents of this paper, it is enough to consider quadratic numbers with purely periodic 

( Si s 2 \ 

we have the equality 

S3 s 4 J 

o S3 + S4H , 

1 1= -— with S1S4 — S2S3 = ±1, 

si + s 2 U 

expressing the equivalence of the number O, with an eventually periodic continued fraction, with the number Q 
with a purely periodic one. Then, it can be shown that our main result (Theorem [T]) applies to both numbers fl 
and for e small enough, and we only need to consider the purely periodic case. For instance, the results for the 
golden number Q = [ 1 ] also apply to the noble numbers Q = [ b ±,..., b r , 1 ]. We point out that the treshold in e 
of validity of the results, not considered in this paper, would depend on the non-periodic part of the continued 
fraction. 

2. This proposition provides a particular case of an algebraic result by Koch |Koc99| . which also applies to higher 
dimensions: for any given vector w 6 R" whose components generate an algebraic number field of degree n , there 
exists a unimodular matrix T having w as an eigenvector with the associated eigenvalue A of modulus greater 
than 1. This result is usually applied in the context of renormalization theory, since the iteration of the matrix T 
provides successive rational approximations to the direction of the vector w (see for instance (Koc99t |Lop02|). 


2.2 Resonant sequences 

In this section and the next one, we review briefly the technique developed in [DG031 for classifying the quasi-resonances 
of a given frequency vector w = (l,f2) whose ratio f l is quadratic, and study their relation with the convergents of 
the continued fraction of Q. A vector k £ Z 2 \ {0} can be considered a quasi-resonance if (fc,w) is small in modulus. 
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To determine the dominant harmonics of the Melnikov potential, we can restrict to quasi-resonant vectors, since the 
effect of vectors far enough from resonances can easily be bounded. 

More precisely, we say that an integer vector k ^ 0 is a quasi-resonance of u> if 

It is clear that any quasi-resonance can be presented in the form 

k°(q) : = (-P, q), with p = p°{q) := rint(gfi) 

(we denote rint (a;) the closest integer to x). Hence, we have the small divisors (k°(q),oo) = qQ — p. We denote by A 
the set of quasi-resonances k°(q) with q > 1 (which can be assumed with no loss of generality). We also say that k°(q) 
is an essential quasi-resonance if it is not a multiple of another integer vector (if p^O, this means that gcd (q,p) = 1), 
and we denote by Ao the set of essential quasi-resonances. 

As said in Section [Til the matrix T = T(f2) given by Proposition [4] (in both cases of purely or non-purely periodic 
continued fractions) provides approximations to the direction of ui = (1,12). Instead of T, we are going to use another 
matrix providing approximations to the orthogonal line (w)- 1 -, i.e. to the quasi-resonances of u. Notice the following 
simple but important equality: 

<(T- 1 ) T fc, w) = (fc, T -1 w) = j <fe,w), (18) 

with A = A(12) as given by Proposition [4] With this in mind, for a quadratic ratio with an (eventually) m-periodic 
continued fraction, we define the matrix 

U = U(fl) := ct(T - 1 ) t , where ct := det T = (-l) m (19) 

(the sign a , which is not relevant, is introduced in order to have a simpler expression in (1^1) 1. It is clear from (fTSl) 
that if k £ A, then also Uk £ A. We say that the vector k = k°(q) = (— p,q ) is primitive if k £ A but U~ 1 k ^ A. 
If so, we also say that q is a primitive integer, and denote V the set of primitive integers, with q > 1. We deduce 
from (1181) that k is primitive if and only if the following fundamental property is fulfilled: 

^<\{kM\<\. (20) 

If a primitive k°(q) = (—p, q) is essential, we also say that q is an essential primitive integer, and we denote Vq C V 
the set of essential primitive integers. 

Now we define, for each primitive vector k°(q), the following resonant sequences of integer vectors: 

s{q,n) := U n k°(q), n > 0. (21) 

It turns out that such resonant sequences cover the whole set of vectors in A, providing a classification for them. 

Remark. A resonant sequence s(q,n) generated by an essential primitive k°(q) cannot be a multiple of another 
resonant sequence. Indeed, in this case we would have k°(q) = cs(q, no) with c > 1 and no > 0, and hence k°(q) would 
not be essential. 


Let us establish a relation between the resonant sequences s(q,n), and the convergents of f2. Alternatively 
to the convergents w(j) = ( qj,Pj) considered in Section [2.11 we rather consider the “resonant convergents” (see 
also |l)Cr;i lb| ). 

v(j) ■■= (~Pj,qj)- 


The next lemma shows that the action of the matrix U defined in dlOl) . on the vectors v(j), is analogous to the action 
of T on the vectors w(j) (which has been described in Proposition U) • This implies that the sequence of resonant 
convergents is divided into m of the resonant sequences defined in (l2lT) . We also see that the primitive vectors 
generating such sequences are the m first resonant convergents (belonging to A). 
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Lemma 5 


(a) Let LI be a quadratic number with an (eventually) periodic continued fraction, LI = [ bi ,..., b r , a\,..., a m ] 
(with r > 0). Then, we have 

Uv(j) = v(j + to), j>r, 

and hence the sequence of resonant convergents v(j), for j > r, is divided into m resonant sequences. 

(b) If LI has a purely periodic continued fraction, Lt = [a\,... , a m ], the primitive vectors among the resonant 
convergents are 

^(1), • • -,v{m) if ai = 1; 

u(0),..., v(m — 1) if a\ > 2. 


Proof. We use the following simple relation between the entries of the matrices T and U (valid in both cases r = 0 
or r > 1): 


if 




then U = 



( 22 ) 


where we have taken into account that detT = (—l) m . Then, the equality Tw(j) = w(j +m), which holds for j > r, 
is exactly the same as Uv(j ) = v(j + m), as stated in (a), by the relation between the vectors v(j ) and w(j). We have, 
as an immediate consequence, that the sequence of resonant convergents v(j) (for j > r) is divided into m resonant 
sequences. 


To prove (b), we first see that the small divisors associated to the resonant convergents v(j) satisfy the equality 

qjLl -Pj = {-l) 3 xo---Xj, j> 0, 

where (xj) is the sequence introduced in (fT71) . This can easily be checked by induction, using the recurrence (fT5l) and 

the equality- aj = Xj . 

x j— i 

If the continued fraction of LI is purely periodic, recalling the expression for A given in Proposition |U(a) and the 
fundamental property (1201) . it is clear that a resonant convergent v (j ) is primitive if and only if the following inequalities 
hold: 

xo * * * Xm—l 1 

---< X 0 ■ ■ ■ Xj < - . 

Recall that Xj £ (0, 1) for any j. Using (a), we see that such inequalities can only be fulfilled by, at most, m consecutive 
values of j. For ai > 2, the first one is j = 0 since Xo = LI < 1/2, and the last one is clearly j = rri — 1. Instead, 
for ai = 1 the first one is j = 1 since xo = LI > 1/2 and XqX\ = 1 — Xq < 1/2, and the last one is j = m since 

Xryi — Xq 1/2. CH 


Remarks. 

1. The matrices T and U cannot be triangular, i.e. we have 6^0 and c 7^ 0 in (1221) . Indeed, this would imply that 
the eigenvalue A is rational, and hence the frequency ratio Ll would also be a rational number. 

2. The primitive resonant convergents given in part (b) of this proposition are all essential primitive vectors, since 
all the convergents Pj/qj are reduced fractions (as a consequence of the fact that the matrices in (fTCl) are 
unimodular). 


2.3 Primary and secondary resonances 

Now, our aim is to study which integer vectors k fit best the Diophantine condition ( 0 . As in mm, we define 
the “numerators” 

7 *:=|<fc,o;)|-|fc|, k£ Z 2 \{0} (23) 
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where we use the norm |-| = |-| L (i.e. the sum of absolute values). As said in Section [2721 we can restrict ourselves to 
vectors k = k°(q) € A (with q > 1 ), and such vectors will be called primary or secondary resonances depending on 
the size of 7 We are also interested in studying the “separation” between both types of resonances. 

Recall that the matrix T given by Proposition 2] has oj = (l,f2) as an eigenvector with eigenvalue A > 1. We 
consider a basis ui, v 2 of eigenvectors of T, where the second vector v 2 has the eigenvalue a/X (of modulus < 1; recall 
that a = detT). For the matrix U defined in (flUl) . let u\, u 2 be a basis of eigenvectors with eigenvalues a/X and A, 
respectively. Writing the entries of the matrices T and U as in (1^1 ) . it is not hard to obtain expressions for such 
eigenvalues and eigenvectors: 


A — a ~\~ bXl, 


?- = d-bSl, 

A 


v 2 = (—bXl,c), ui = (c,bXl), zt 2 = (—1). 

We also get the quadratic equations for the frequency ratio and the eigenvalue A: 

bXl 2 = c — (a — d) fl, A 2 = (a + d) X — a. 

For any primitive integer q , recalling that we write k°(q) = (—p,q), we define the quantities 

r q := ( k°(q ), u) = qXl - p, z q := (k°(q), v 2 ) = cq + bpXl. 


(24) 


(25) 


(26) 


Remark. As a consequence of the fact that XI is an irrational number, one readily sees that, if q ^ q, then r q ^ r q 
and z q ^ z q (in the latter case, using also that c / 0, as seen in remark [T| after Lemma [5]). 


The following proposition, whose proof is given in l)( 11)3 (see also IDGG14a1 for a comparison with the case of 
3-dimensional cubic frequencies), shows that the resonant sequences s(q,n ) defined in (1211) have a limit behavior: the 
sizes of the vectors s(q,n) exhibit a geometric growth , and the numerators 7 s ( q ,n) tend to a “limit numerator” 7 *, 
as n —► 00 . 


Proposition 6 Let uj = (1, XT) be a frequency vector whose ratio £ (0,1) is quadratic, and consider the vectors v 2 
and u 2 as in f 2/\ ). For any primitive integer q £ V (see 12(H) ). consider the quantities r q and z q defined in 126]) . Then, 
the resonant sequence s(q,n) defined in (2l\) satisfies: 


(a) \s(q,n)\ = K q X n + 0(X n ), where K q := 


(u 2 ,v 2 ) 


u 2 


(b) 7 s(q,n) =1 q +0(X 2n ), where 7 * := lim 7 s(<? , n) = \r q \ K q . 


Using (124H261) . we get the following alternative expression for the limit numerators: 

7 « = I 59 ' ’ wher e 6 q :=^L = cq 2 - (a-d)qp-bp 2 , p=p°(q). (27) 

It is clear that S q ^ 0 and it is an integer. We can select the minimal of the values \5 q \ and, consequently, of the limit 
numerators 7 *, which is reached by some concrete primitive q. We define 

S * '' = fev ^ = 1^1 - 1 ' 7 * := = b'l > °- ( 28 ) 

It is easy to see, as a consequence, that lim inf 7 *, = 7 * > 0. Hence, any vector with quadratic ratio satisfies the 

|fc|—»oo 

Diophantine condition (f7|) . and we can consider 7 * as the asymptotic Diophantine constant. 

As we see, all limit numerators 7 * are multiple of a concrete positive number. An important consequence of this 
fact is that it allows us to establish a classification of the vectors in A. We define the primary resonances as the 
integer vectors belonging to the sequence So(n) := s(q,n ), and secondary resonances the vectors belonging to any of 
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the remaining sequences s(q,n ), qf=q (recall that q is the primitive giving the minimum in ([28]) ]). We also introduce 
normalized numerators jk and their limits 7 *, q £ V, after dividing by 7 *, and in this way 7 / = 1 . We also define a 
value B 0 = B o(I2) measuring the “separation” between the primary and the essential secondary resonances: 


7 k ■= 



-* . = t£ = N 

lq ' 7 * <5* 5 


B 0 := 


min 

q£Vo\{q} 


iq- 


(29) 


Using the fundamental property (l20l) and the inequality \p — gf2| < 1/2, we get the following lower bound for the 
limit numerators, which sligthly improves the analogous bound given in |DG03] : 


7? > 


(1 + f l)q - a 

2A 


|6|U(1 + D) 

2 |c + bfl 2 \ 


(30) 


Remarks. 

1. Since the lower bound (13U1) is increasing with respect to q 1 it is enough to check a finite number of cases in order 
to find the minimum in (1281) . 

2. We are implicitly assuming that the primitive integer q providing the minimum in (1281) is unique. In fact, we 
will show in Section 12.41 that this is true for the cases of metallic or metallic-colored ratios ft, introduced in 
Section 11.31 But in other cases, the minimum could be reached by two or more primitives and, consequently, 
there could be two or more sequences of primary resonances. For instance, it is not hard to check that for the 
ratio Q = [ 1 , 2 , 2 ] there are two sequences of primary resonances. 

3. Any primitive integer q generating a sequence of primary resonances is essential. Indeed, if q is not essential, 
then we have k°(q) = cs(q,no) with c > 1 and no > 0 , and therefore s(q,n ) = cs(q,no + n), which implies 
by (l23l) that = c 2 y|, and the minimum in (l28l) would not be reached for q. 


Next, we show that the sequence of primary resonances is one (or more) of the m resonant sequences in which, by 
Lemma 0 the resonant convergents are divided if the continued fraction of f 1 is m-periodic. In fact, we can give a 
lower bound for the numerators of all the remaining sequences. 


Lemma 7 For any primitive integer q such that the vectors in the sequence s(q,n ) are not resonant convergents, its 
normalized numerator satisfies 7 * > \Jhj1. 


Proof. We use some results in |Sch801 §1.5] (namely, Theorems I.5B and I.5C), concerning the properties of the 
convergents of any irrational number. On one hand, for an infinite number of convergents the inequality \q n Q — p n \ < 
1/(s/5 q n ) is satisfied; and on the other hand, if a given integer q > 1 is not a convergent and p/q is a reduced fraction, 
then |gfi — p\ > 1/2 q. To compare such results with our Diophantine condition J7|, notice that |fc°(q)| = q+p « (l+fi )<7 
as q —> 00 . 

The first quoted result implies that, at least for one of the resonant sequences s(q,n) whose vectors are resonant 
convergents, its limit numerator satisfies 7 * < (1 + S2)/-\/5. By the second result, if a given resonant sequence s(q,n ) 
is generated by an essential primitive q and its vectors are not resonant convergents, then 7 * > (1 + S2)/2. This is also 
true if q is not essential, by the previous remark 0 Dividing the two bounds obtained, we get the lower bound -s/5/2 
for the normalized limit 7 *, when q does not generate a sequence of resonant convergents. □ 


2.4 Results for metallic and metallic-colored ratios 

Now, we provide particular arithmetical results for the (purely periodic) cases of a metallic ratio fl = [a], and a 
metallic-colored ratio D = [ a, b ], introduced in Section 11.31 
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Metallic ratios. Let us write, for a given f 1 = [a], a > 1, the matrix T = T(D) and the eigenvalue A = A(f l), as 
deduced from Proposition 0|a), and the matrix U = t/(D) from (I221l . 


T = 


a 

1 




We also have from (l25|) the quadratic equation 


A 2 = aA + 1. 


A = n=“ 


d. 


(31) 


(32) 


By Lemma O all resonant convergents v(j) belong to a unique resonant sequence, whose primitive vector is v(l) = 
(— 1 ,1) if a = 1, and u(0) = (0,1) if a > 2. We deduce from Lemma[T]that this resonant sequence provides the primary 
resonances: in both cases q = 1 and hence so(n) = s(l,n). In the next result, we compute the separation Bo, defined 
in (|29[) . for all metallic ratios, providing in this way a sharp lower bound for the normalized numerators of all the 
essential secondary resonances. 


Proposition 8 Let D = [a], a > 1, be a metallic ratio. Then, the sequence of primary resonances is generated by the 
primitive integer q = 1, and we have: 


Bo = 7 


* 

Qi 


5 if a = 1, 
a if a> 2 , 


for qi 


7 if a = 1, 

3 ifa = 2, 

ail if a> 3. 


Proof. We use the expression (1271) . taking into account the entries of the matrix T given in (TSTT) . For the primary 
resonances, we have 5\ = 5* = 1, and hence 7 * = • Dividing (13U1) by 7 * and using that A = 1/D we get, for 

the normalized numerators 7 * = |<5 g |, the following lower bound: 


_* 1 + D 2 

7, > 


d 

I ' 


If a = 1 (the golden ratio), one checks that the second essential primitive is (—4, 7) with 77 = 5, and 7* > 5 for q > 8 . 
For a > 2, assuming that q > 2/D we get 7* > a. Otherwise, if q < 2/D, since p < gD — 1/2 we get p = p°(g) < 3/2, 
i.e. p = 0 or p = 1. The only essential primitive with p = 0 is (0, 1), which gives the primary resonances, and for p = 1 

we have an “interval” of primitives (— 1, g), with —i— < g < — and g ^ a (we have applied (1201) together with the 


fact that a < 1/D < a+1). For such primitives, applying G71) we obtain 5 q = q 2 — ag— 1, a quadratic polynomial in g, 
which is a increasing function for q > a/2, with S a ± 1 = ±a. This change of sign indicates that 7 * = |<5 g | is minimal 
for q = a ± 1. This argument is valid for a = 2 (the silver ratio), but in this case we must exclude q = a — 1, which 
lies outside the interval considered. □ 


Metallic-colored ratios. Now, we consider D = [ a, b], 1 < a < b. Recall that, for a = 1, this is called a 
golden-colored ratio; we see below that our results are somewhat different for this particular case. We have 


T = 


ab+1 

b 



( 1 -M, 

y —a ab -\-1 J 


A — --— — ab ~b 1 H - cl£1 

1 — ail 


and, from (f25|) , the quadratic equation 

A 2 = (ab + 2)A - 1. (33) 

Applying Lemma [E] we see that the resonant convergents v(j) are divided into 2 resonant sequences, whose respective 
primitive vectors are 

«(!) = (-!>!)> v(2) = (~b,b + l) if a = 1 ; 

(34 

u( 0 ) = ( 0 ,l), u(l) = (—l,a) if a > 2 . 

By Lemma [3 one of the 2 sequences of resonant convergents provides the primary resonances: so{n) = s(q,n). We 
call the main secondary resonances the vectors in the second sequence, which we denote as Si(n) := s(g, n). Next, 
we find in Proposition [9] the value of yjt, i.e. the limit numerator of the main secondary resonances, and this gives 
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a rigorous upper bound for the separation Bq. In Numerical Result [TO] we present numerical evidence of the exact 
equality Bq = yi, after checking it for a large number of ratios. We point out that, for a given concrete frequency 
ratio 12, the separation Bq(LI) can be rigorously determined since, by the lower bound (1301) . it is enough to consider 
the limit numerators 7 * for a finite number of essential primitive integers < 7 . 


We also find in Numerical Result m the value of the “second separation” , i.e. the minimal normalized limit 
numerator among the essential resonant sequences whose vectors are not resonant convergents: 


B 1 


min 

qeVo\{q,q} 


iq- 


(35) 


Proposition 9 Let 12 = [a, b], 1 < a < b, be a metallic-colored ratio. Then, the sequences of primary resonances and 
main secondary resonances are generated, respectively, by primitive integers q, <7 given by 

<7 = 1, g = 6 + l if a = 1 , 
q = a, <7 = 1 if a > 2 . 

In both cases, the separation satisfies 

1 < B 0 < % = - . 

H a 

Proof. We use (1271) in order to determine which primitives (1341) generate the sequence of primary resonances. For 
a = 1, we obtain <5i = —1 and <5{,+i = 6 . For a > 2, we obtain 61 = b and S a = —a. In both cases, the minimum (in 
modulus) is <5* = a, which is reached for <7 = 1 if a = 1, and <7 = a if a > 2. Then, we have <7 = 6 + 1 ifa = 1, and <7 = 1 
if a > 2, and we obtain y£ = \bq/5f\ = 6 /a, which provides the upper bound: Bq < 6 /a. We also have the inequality 
B 0 > 1 as a consequence of Lemma [3 □ 


Numerical Result 10 In the notations of Proposition^ the separation and the second separation are exactly 0 

6 + 4 if a = 1 , 


Bo = a =- 

y a 


B 1 = 


(a — 1)6 + 1 


if a > 2 , 


(36) 


which satisfy 1 < Bq < Bi. 


Justification. Numerically, we can compute B\ by bounding from below the limit numerators 7* for all the essential 
primitives <7 < 7 , <7 (in view of (l30l) . only a finite number of primitives q have to be considered). We have checked that 

they all satisfy 7* > 6 /a (at least for all the frequency ratios we have explored), and hence we get B 0 = b/a , and 
B\ > B 0 . We also obtain an expression for B 1 , given in (136)) separately for the cases a = 1 and a > 2. □ 

Remark. The numerical explorations allow us to determine accurately the primitive integers <72 such that B\ = y * 2 , 
i.e. giving the minimum in (1351) : 


2 

if a = 1, 6 = 2 , 

3, 9 

if a = 1 ,6 = 3, 

4, 11 

if a = 1, 6 = 4, 

5, 8 , 13 

if a = 1, 6 = 5, 

6+3 

if a = 1,6 > 6 , 

26+3 

if a = 2 , 

a— 1, a 6 + a+ l 

if a > 3. 


In each case, the primitive integers <72 generate the “third most resonant” sequences among the non-convergent ones 
(i.e. after the 2 sequences of resonant convergents). Again, we stress that it is possible to obtain this kind of results 
thanks to the lower bound (1301) . which allows us to carry out a finite number of computations for any given ratio fl. 


3 The values of Bq and B\ have been checked numerically for all golden-colored ratios with 1 = a < b < 10 6 , and for all metallic-colored 
ratios with 2 < a < b < 10 3 . 
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3 Searching for the asymptotic estimates 


In order to provide asymptotic estimates for the splitting, we start with the first order approximation, given by 
the Poincare-Melnikov method. Although our main result (Theorem [L} is stated in terms of the splitting function 
A4(6) = V£(0), it is more convenient for us to work with the (scalar) splitting potential C(8), whose first order 
approximation is given by the Melnikov potential L(0). 

In this section, we provide the constructive part of the proof, which amounts to find, for every sufficiently small e, 
the first and the second dominant harmonics of the Fourier expansion of the Melnikov potential L(9 ), with exponentally 
small asymptotic estimates for their size, given by functions hi(e) and / 12 (e) in the exponents. As a direct consequence 
of the arithmetic properties of quadratic ratios, such functions are periodic with respect to lne. We also study, 
from such arithmetic properties, whether the dominant harmonics are given by primary resonances. This allows us 
to provide a more complete description of the functions hi (e) and h? (e) in some particular cases (Theorem [2] and 
Numerical Result 0 . 

The final step in the proof of our main result is considered in Section [I] It requires to provide bounds for the 
sum of the remaining terms of the Fourier expansion of L(9 ), ensuring that it can be approximated by its dominant 
harmonics. Furthermore, to ensure that the Poincare-Melnikov method (j2j) predicts correctly the size of the splitting in 
the singular case /1 = e r , one has to extend the results to the Melnikov function M{9) by showing that the asymptotic 
estimates of the dominant harmonics are large enough to overcome the harmonics of the error term in (0. This step 
is analogous to the one done in |DG04I for the case of the golden number f2 = [ 1 ] (using the upper bounds for the 
error term provided in |DGS04n . 


3.1 Estimates of the harmonics of the splitting potential 


We plug our functions / and h, defined in into the integral 
potential, where the coefficients can be obtained using residues: 

L{6) = Lk cos (( fc ’ °) ~ CT fc), 

fc€2\{0} 

We point out that the phases 07 are the same as in ©• Using m and taking into account the definition of the 
numerators 7 k in (1231) . we can present each coefficient Lk = Lk(s), k £ Z \ {0}, in the form 

L* = «e-' 5 -, a * (£)cs |yl- * (£) = ' ,W + 27775’ (37) 

where an exponentially small term has been neglected in the denominator of 07 . The most relevant term in this 
expression is 0k, which gives the exponential smallness in e of each coefficient, and we will show that 07 provides a 
polynomial factor. This says that, for any given e, the smallest exponents 0k(s) provide the largest (exponentially 
small) coefficients Lk{e) and hence the dominant harmonics. We are going to study the dependence on e of this 
dominance. 


(fl3|) and get the Fourier expansion of the Melnikov 


_ 27r|(/c,a; e )| e~ p ' fc ' 
sinh \ £(k.ujr)\ 


We start with providing a more convenient expression for the exponents 0k(z), which shows that the smallest ones 
are 0(£ -1 / 4 ) (this is directly related to the exponents 1/4 in Theorem [T|. We introduce for any given X, Y the 
function 


G(e-,x,Y) 



(38) 


which has its minimum at e = X with G(X; X, Y ) = F 1 / 2 as the minimum value. Notice that each function G(-; X, Y) 
is determined by the point (X, Y 1 / 2 ). Now, we define 


9 k(e) ■■= G(e;efc, 7 fc ), 


£fc := 


D o lk 

| fc | 4 


D 0 := 



19 













and the functions < 7 fc(e) have their minimum at e = £&, with the minimal values gk{^k) = 7 %. Recall that the 
asymptotic Diophantine constant 7 * = 7 ~ and the normalized numerators 7 & = 7 ^/ 7 * were introduced in (128112911 . We 
deduce from (l37l) that 


Pk(e) = ~Tfl Co := (2ttp7*) 1/2 , 


(39) 


and hence the lower bound 


Pk{e) > 


cwT! 

e 1 / 4 


Since we are interested in obtaining asymptotic estimates for the splitting and its transversality, rather than lower 
bounds, we need to determine for any given e the first and the second essential dominant harmonics, which can be 
found among the smallest values <?&(£). To this aim it is useful to consider, for a given frequency ratio f2, the graphs of 
the functions gk{s) associated to essential quasi-resonances k £ Aq (recall that the notion of “essentiality” has been 
introduced at the beginning of Section [272]) . As an illustration, such graphs are shown in Figure Oja)] for a concrete 
example (the bronze ratio = [3]), using a logarithmic scale for e. Other examples are shown in Figures [2][3] 
The periodicity which can be noticed from the graphs can easily be explained from the classification of the integer 
vectors into resonant sequences (recall their definition in (1211) 1. Indeed, for k = s{q,n) belonging to a concrete 
resonant sequence, using the approximations for \s(q, n)| and 7 s ( q , n ) given by Proposition [ 6 ] we obtain the following 
approximations as n — > 00, 


gs(q,n){e) ~ g*s(q, n)( £ ) G{ £ \ £ l(q,n) ) 7q), 


£ s(q,n) ~ £ s(q,n) 


^(T'g) 2 

AT 4 A 4 " 


(40) 


which motivates the use of a logarithmic scale. We point out that the graphs shown in Figure [ f)7i)| do not correspond to 
the true functions g s (q, n )( £ )i but rather to the approximations „)(£), which satisfy the following scaling property: 


9* s ( q ,n+i)( £ ) =9*s( q ,n)(^z)- ( 41 ) 

This gives, for any resonant sequence, the mentioned periodicity: the graph of nJt i\ is a translation of g*^ q n y to 
distance 4 In A. For non-essential resonant sequences, whose vectors do not belong to Ao, we see that, if s{q,n) = 
cs(q , no + n) with c > 1 and no > 0 , then 


9*s(q,n)( £ ) — C 9*s(q,n 0 +n)( £ ) 

(see also the remark [3] just before Lemma [7}. 


(42) 


In order to study the dependence on e of the most dominant harmonics, it is useful to study the intersections 
between the graphs of different functions g%.{s), since this gives the values of £ at which a change in the dominance 
may take place. In the next lemma, we consider the graphs associated to two different quasi-resonances k,k £ A, and 
we show that only two situations are posible: they do not intersect (which says that one of them always dominates 
the other one), or they intersect transversely in a unique point (and in this case a unique change in the dominance 
takes place). 


Lemma 11 Let k,k £ A, with k ^ k, given by k = s{q,n) and k = s(q,n), and assume that 7* < 7^ . Denoting 
Z = [eL/e^ and W = ( 7 '|/ 7 g) 1 ^ 2 ) the graphs of the functions gl(e) and g^e) intersect if and only if Z < 1/W 

or Z > W. If so, the intersection is unique and transverse, and takes place at e = £^ • 


/Z(WZ-1)\ 2 

V Z-W J ‘ 


Proof. First of all, we show that gl and g^~ cannot be the same function. By the definition (1551) . if g%. = g^~ then 
we have 7 * = 7 ^ and £^ = £ j- The latter equality implies that K q X n = Kq A”. Using the expressions given in 
Proposition [Gl we get the equalities \r q z q \ = \r q z q \ and z q \ n = z q A”. We deduce that \r q /r q \ = A n ~ n , but we have 
\r q \ , \r q \ £ (1/2A, 1/2) by the fundamental property (l20l) . This says that n = n and hence z q = z q . As seen in the 
remark next to the definition (12611 . we also get q = q, which contradicts the assumption k k. 


Now, introducing the variable = (e/e ^) 1 ^ 4 > 0, we define 


/i(C) := 


9k( £ ) 

V (%*) 1/2 


C+7 = 


/ 2 (C) := 


w ( c z 
~2 [z + C 


4( £ ) 

(%*) 1/2 
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with W > 1 by hypothesis, and it is clear from the above analysis that we cannot have W = Z = 1. It is straightforward 

to check that the graphs of f± and /2 can intersect only once, transversely, at ( 2 = —--—— . Such an intersection 

Z — vv 

occurs if and only if Z < 1/W or Z > W. Then, we get the result after translating from £ to the original variable e. 

□ 


The sequence of primary resonances so(n) = s(q,n), introduced in Section 12.31 plavs an important role, since 
they give the smallest minimum values among the functions ^(e), and hence they will provide the most dominant 


harmonics, at least for e close to such minima. With this fact in mind, and recalling that = 1, we denote 


3n(e) : = 9*s 0 (n)( £ ) = G(e;£„, 1) = - 


1/4 


£n •— £ 


s 0 (n) 


Do 

K~\ Ar 


(tr 


(43) 


To study the periodicity with respect to lne, we introduce intervals I n whose “length” (in the logarithmic scale) 
is 4 In A, centered at e n , and the left and right “halves” of such intervals, 


T •= 

- L ' n • 


'n+ 1 ? £ 


t] -^n ? -^n • [^n+l^n] ? n * ' 


(44) 


where e' n := yje n e n ^\ = \ 2 e n are the geometric means of the sequence e„. For a given n > 1, it is easy to determine 
the behavior of the functions (143)) : for e G I+, the value of the function g n (e ) decreases from Ji to 1, the value of the 
function ~g n+ i(s ) increases from J\ to J 2 , and we have g m (e) > J2 if m ^ n,n + 1 (recall that the values Ji and J 2 
were defined in (fl4)l l. A symmetric result holds for e € I~ with the functions g n (e) and g n _ x (e) (see the red graphs 
in Figure QJl'a)| for an illustration). 


3.2 Dominant harmonics of the splitting potential 

In this section, we introduce the functions hi(e) and /12(e) appearing in the exponents in Theorem [lj as the first and 
second minima, for any given e, of the values g^{s) among the essential quasi-resonances k G Aq. We study some of the 
properties of h\(e) and /12(e), which hold for an arbitrary quadratic ratio fl. In Section 13.31 we put emphasis on the 
dependence of such functions on the continued fraction of the frequency ratio f 2 , giving a more accurate description of 
them, for the cases of metallic and metallic-colored ratios, whose arithmetic properties have been studied in Section l^Tfl 

Namely, we provide information on the minimum and maximum values of the functions h\(e) and /12(e), and show 
that they are piecewise-smooth and 4 In A-periodic in lne, and give lower bounds for the number of their corners 
(i.e. jump discontinuities of the derivative) in any period, say the interval l n = [e^.^e^J, or rather the semi¬ 
open interval (e^.^eJJ, to avoid repetitions if the endpoints are corners. In fact, such properties are clear from 
Figures EEH3 for the concrete frequency rations considered there, but we are going to show that they hold for an 
arbitrary quadratic ratio Cl. 

Previously to this, let us define two functions analogous to h\{e) and /12(e), but taking into account only the 
primary resonances: 

Me) : = min 5„(e) = g* Nl {£), Me) : = min 9* n {e) = g* N2 (e), (45) 

n 1 n^N 1 

with Ni = Ni(e). In other words, the two dominant harmonics among the primary resonances correspond to 

~Si = )5j(e) = so(Ni), i = 1,2. 

On each concrete interval I n (see the definition (1441) 1 one readily sees, from the properties described in the last part of 
Section [3.11 what primary resonances provide the first and second minima: iVi(e) = n for e G X n , and N-i (e) = n± 1 
for e G Z„ . It is also clear that the functions /11(e) and /12(e) are piecewise-smooth and 4 ln A-periodic. In each period, 
the function h\{e) has exactly 1 corner (at e' n ), and /12(e) has exactly 2 corners (at e' n and e„). Moreover, we have 

min/ii(e) = /ii(e„) = 1 , max/i 2 (e) = /i2(e„) = J2, 
max/ii(e) = min /12(e) = hi(e' n ) = h 2 (e' n ) = Ji 
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(see also Figure t jfb)! for an illustration). 

Now, we define the functions hi(e) as the minimal values of the functions g k (s) among all essential quasi-resonances, 
and we denote S', = Si(e) the integer vectors k at which such minima are reached: 

hi{e) := min g* k (e) = g* Sl (e), h 2 (e) := min g%(e) = g* s (e), 

keA 0 fce-A 0 \{Si} 

We) := KJ „Tl s n ,A) 9j( ' ) = * (e) ' 

It is clear that hi{e) < hi(e) for any e and i = 1,2. In order to provide an accurate description of the splitting and 
its transversality, we have to study whether the equality between the above functions can be established for any value 
of e, or at least for some intervals of e. This amounts to study whether the dominant harmonics can be always found 
among the primary resonances (Si = Si) or, on the contrary, secondary resonances have to be taken into account. 

In fact, the properties described above for the functions /q(e) are partially generalized to the functions hi(e) in the 
next proposition, which corresponds to some parts of the statement of Theorem [TJ concerning such functions. Recall 
that the values J\ and J 2 were defined in (fl4l) . 


Proposition 12 The functions hi(e) and h 2 (e) are piecewise-smooth, 41n X-periodic in lne. In each period, the 
function h\(e) has at least 1 corner and h 2 (e) has at least 2 corners. They satisfy for e > 0 the following bounds: 

min h\ (e) = 1, max hi(e) < J\, max h 2 (e) < J 2 , hi(e) < h 2 (e). 

The corners of h\(e) are exactly the points e such that h\(s) = h 2 (e). The corners of h 2 (e) are the same points e, and 
the points e where / 12 (e) = / 13 (e). 


Proof. First of all, it is clear that the functions h\ and h 2 are 4 In A-periodic, as we see from the scaling property (HQ- 
Then, we can restrict ourselves to a concrete interval, say I±. Recalling also that /q(e) < /q(e) < Ji for any e and 
i = 1 , 2 , the minimum in the definition (14611 of hi can be restricted to the integer vectors k = s(q,n) £ Ao such that 
the graph of the function g*^ q n ^ visits (the interior of) the rectangle T\ x [ 1 , Jf. We are going to show that this is 
possible only for a finite number of integer vectors. 

Indeed, recalling (l38l) . if the graph of a function G(e\X, Y) visits I\ x [1, Ji] then F 1 / 2 < J x and e 2 < X < £ 0 ; 
and if it visits I\ x [1, J 2 ] then F 1 / 2 < J 2 and £3 < X < e' 0 . For the function g* (q ,, defined in (|40l) . we have to 
consider F = 7 * and X = e*^ q . By the lower bound (Idlil) , it is clear that only a finite number of functions g*^ q n j 
can visit the rectangle Ii x [1, Jf[, i = 1 , 2 . This implies, by Lemma fTTl that only a finite number of (transverse) 
intersections between the graphs of g*, n ^ can take place inside the rectangles Ii x [ 1 , 

We deduce from the above considerations that the functions h \ and h 2 are piecewise-smooth. Indeed, we can 
consider a partition of T\ into subintervals such that, for £ belonging to (the interior of) each subinterval, the function h\ 
coincides with only one of the functions g*^ q n y i.e. the dominant harmonic is given by S’i(e) = s(q,n), which remains 
constant on this subinterval. At each endpoint of such subintervals, a change in the dominant harmonic takes place, 
i.e. Si has a jump discontinuity. By Lemma fill the endpoints of the subintervals correspond to transverse intersections 
between the graphs of different functions g*^ q n y which give rise to corners £ of hi. A similar argument applies to the 
function h 2 , with a different partition, associated to the changes of the second dominant harmonic S 2 (e). In fact, the 
values £ are the points where hi(e) = h 2 (i ), and they are corners of both functions hi and h 2 . In the same way, the 
function h 2 has additional corners e at the points where h 2 (e) = / 13 (E). 

Finally, we provide a lower bound for the number of corners £, e in a given period (if the endpoints of a period are 
corners, we count them as one single corner). Since 3 i(£i) = 1, we have S\(e) = S'i(e) = so(l) in some neighborhood 
of £1 £ T\. Analogously, we have S\(e) = 51 (e) = so( 2 ) in some neighborhood of e 2 £ I 2 , which implies the existence 
of at least one corner of hi with e 2 < e < £1 and, consequently, in any given period. On the other hand, if £ < £ 
are two consecutive corners of hi (and h 2 ), there exists at least one additional corner e of h 2 (and / 13 ), since S i(e) 
and 5 2 (e) cannot be simultaneously constant in the interval [e, e] (this would imply that gg and </£ intersect at both 
points £ and £, which is not possible by Lemma fTTl) . □ 

Remarks. 
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1. We can also deduce from the proof of this proposition some useful properties of the functions Si = Si(e), giving 
the dominant harmonics. Namely, each function 5 j(e) is “piecewise-constant”, with jump discontinuities at the 
corners of hi(e). Moreover, the asymptotic behavior of the functions Si(e) as e —> 0 turns out to be polynomial: 

ISOOI-^TI- ( 47 ) 

Indeed, the first dominant harmonic belongs to some resonant sequence: we can write S'i(e) = s(q, N) for some 
q = q(e), and for N = N(e) such that the value e*^ q N -> is the closest to e, among the sequence e*, q n y n > 0. 
Recalling (HOI) and the estimate |s(g, N)\ ~ X N given in Proposition [SI a), we get (H71) . An analogous argument 
holds for 52 (e), possibly replacing N by N ±1, and possibly belonging to a different resonant sequence s(q, •). 
Notice that it is not necessary to include q in the estimate (UTl) (in spite of the fact that K q and 7 * appear in the 
expression OT), since by the arguments in the above proof (Proposition 1121) only a finite number of resonant 
sequences s(q , •) can be involved. 

2. A more careful look at the arguments of the previous remark, says that, if the dominant harmonics in a given 
interval X n are known, then in view of the scaling property m the dominant harmonics in the interval I n +i 
are the next vectors in the respective resonant sequences: 

Si(e) = USi( A 4 £ ) 

(recall that the matrix U appears in the definition of the resonant sequences in (1211) 1. 

3. Although we implicitly assume that there exists only one sequence of primary resonances (see remark [2] before 
Lemma Q, it is not hard to adapt the definitions and results to the case of two or more sequences of primary 
resonances. In this case, we would choose one of such sequences as “the” sequence So(n), when the functions 
g n (e) are defined in (ld3l) . As an example, we show in Figure [3] the graphs of hi(e) and /12(e) for the ratio 
Cl = [ 1, 2, 2], with two sequences of primary resonances. 


3.3 Dominant harmonics for metallic and metallic-colored ratios 

This section is devoted to the proof of Theorem [2] and the justification of Numerical Result [3j providing a more 
accurate description of the functions h\{e) and / 12 (e) for the cases of 1-periodic and 2-periodic continued fractions, 
i.e. for metallic and metallic-colored ratios Cl, introduced in Section [ 1.31 using some arithmetic results from Section fe.dl 
We emphasize the different behavior of the two functions in each case. 




Figure 3: Graphs of the functions g*( 9rt )(e), /ii(e) and / 12 (e) for Cl = [1,2,2]: the two sequences of primary resonances 
correspond to the red and magenta graphs; the non-essential resonances are not represented. 
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The main issue is to discuss whether the first dominant harmonic S'i(e) and (eventually) the second one S 2 (e) are 
given, for any e, by primary resonances: S'i(e) = Si(e). If so, the function /ii(e) and (eventually) the function / 12 (e) 
coincide with the functions hi(e) introduced in (1451) . whose description is very simple (as in Figure [[]). Otherwise, 
the dominant harmonics are given by secondary resonances at least for some intervals of e, which leads to a more 
complicated function / 12 (e) (than / 12 (e), as in Figure | ^a)[ ), or both complicated functions /ii(e) and / 12 (e) (as in 
Figures QKb)| and [3]). 

The proof of such results requires a careful analysis of the role of secondary resonances, is carried out rigorously for 
the statements of Theorem [2] Instead, for the statements of Numerical Result [3] we provide an evidence after having 
checked them numerically for a large number of frequency ratios. 

Proof of Theorem\2\(a) / Justification of Numerical Result\f^a) . Both of these results concern the behavior 
of hi(e) for a metallic ratio or a golden-colored ratio. In the first case we provide a rigorous proof, and in the second 
case the result relies on Numerical Result m which has been validated numerically for a large number of cases. 

It will be enough to show the lower bound 

> Ji, (48) 

where Bq is the separation between the primary and the essential secondary resonances (recall the definition (1291) 1. 
This lower bound ensures that the most dominant harmonic is found, for all e, among the primary resonances: 
S i(e) = S i(e), and hence h\(e) = h\{e) (such facts are reflected in Figures [ J][b)| and E jla)[ ), which completes the proof, 
in view of the properties of the function hi defined in (1451) . 

Hence, it remains to show that the inequality (l48l) is fulfilled in the two cases of a metallic and a golden-colored 
ratio. Notice that, by the definition of Ji in (TTH) . we can rewrite C51) as 4B 0 A > (A + l) 2 . 

For a metallic ratio 17 = [a], a > 1, we know from Proposition [5] that B 0 = 5 if a = 1, and B 0 = a if a > 2 
(a rigorous result). Then, the inequality C51) is easily checked using the quadratic equation (1351) . 

On the other hand, for a golden-colored ratio 17 = [1,6], b > 2, by Numerical Result [TOl we have Bq = b (which 
has been established for 2 < 6 < 10 6 ). Then, it is easy to check the inequality (1481) using in this case the quadratic 
equation (l33l) . □ 

Proof of Theorem\^(b ). Let us consider a golden-colored ratio 17 = [ 1, 6], b > 2. We know from Proposition [9] 
that the primary resonances so(n) = s(q,n ) and the main secondary resonances s\(n) = s(q,n) are generated, respec¬ 
tively, by q = 1 and q = b + 1, or, equivalently, by the vectors u(l) = (—1,1) and v{2) = (—6, b + 1). To study the 
relative position of the graphs of the functions 3* l( - ra ) with respect to the functions ~g n = «/*,*, we compute: 

g : i( n) (7 b * +1 fl<t 6 2 (A — l) 4 1 

£n K± +l (6A) 4 A 2 ’ 

where we have used the quadratic equation (1531) and the fact that, by (1551) . 

-Kfc- )-i Zb -|-i 6(6 + 1) + 617 6A 

Ki ~ — ~ 6 + 17 “ A- 1 ' 

Hence, we have seen that = e n /A 2 = s' n _i, he. the geometric means introduced in (TH1) (see also Figure [ Ja)] ). 

Now, let us check that 

\JBq Ji < Ji- 

We know from Proposition [9] that B 0 < 6 (a rigorous result). Then, it is enough to see that 
which can be easily checked using again the quadratic equation (l33l) . 

Notice that g Sl (n)(Sn) = 3si(n+i)(£n) = \fBoJi- We deduce from (l49|) that, for some interval around z n . the 
second dominant harmonic Si is not the primary resonance Si = so(n — 1) (if e > e n ) or Si = sg(n + 1) (if e < e n ), 
since at least a main secondary resonance is more dominant: si(n) (if e > e n ) or si(n + 1) (if e < e n ) (eventually, 
another secondary resonance could also be the second dominant harmonic Si). This implies that at least 3 changes in 


(49) 

6 A(A + 1) 2 < (A 2 + l) 2 , 
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the second dominant harmonic take place in a given period, and hence the function /12 has at least 3 corners. We also 
deduce that the maximum value of the function /12 is < J 2 , since the value J 2 can only be reached, at the points e n , 
if a primary resonance is the second dominant there (again, see Figure QK a) | for an illustration, where hi has 4 corners 
in each period). 

Concerning the minimum of the function h 2 , it is always reached at the points by a primary resonance: 
g n (e' n ) = 9 n-i(^n) = since by the inequality (1551) all secondary resonances take greater values at e' n . □ 

Proof of Theorem\^(c). Now, we consider a metallic-colored but not golden-colored ratio, f2 = [a, 6 ], 2 < a < b. 
In this case, we know from Proposition [9] that the primary and main secondary resonances, so(n) and si(n), are 
generated, respectively, by q = a and q = 1, or, equivalently, by the vectors u(l) = (—l,a) and u(0) = (0,1). As in 
part (c), we study the relative position of the functions g* Sl ^ with respect to g n , by computing: 

<(n) = mfKg 2 = (V«) 2 (A-1) 4 = 2 

£n K? b 4 

where we have used the quadratic equation (1331) and the fact that 

K a z a ab + afl A — 1 
~K X = ~z x = b = b ' 

Hence, we have seen that £*^ n \ = \ 2 e n = e' n (see also Figure [ jfb)] ). 

Next, we show that, instead of (1451) . we have 

\Z~Bq < J\ 

or, equivalently, 4U 0 A < (A+ l) 2 . We know from Proposition [9] that B 0 < b/a (a rigorous result). Then, it is enough 
to see that 4b A < a(A + l) 2 , which can be checked using again the quadratic equation (1351) . 

We deduce that, for some interval around e ' n , the most dominant harmonic is not a primary resonance: Si / ^i, 
since at least the secondary resonance si(n) is more dominant. This implies that at least 2 changes in the dominance 
take place in a given period, and hence the function h\ has at least 2 corners. We also deduce that the maximum value 
of the function hi is < Ji, since the value Ji can only be reached at the points e ' n , provided a primary resonance is 
the most dominant there (again, see Figure [^b)] for an illustration). □ 

Justification of Numerical Result{3\(b) . We consider a metallic ratio V. = [a], a > 1, and we are going to show 
that, for any e, the second dominant harmonic is also a primary resonance: Si(e) = 82 (e), and hence /j- 2 (s') = hi(e) 
(see Figure [T] for an illustration). Then, it is enough to use the simple properties of the function /12 defined in (|45jl . 

Thus, we have to check that a secondary resonance cannot be the second dominant harmonic in any interval 
of e. By the periodicity, we can restrict ourselves to primitive vectors: s(q, 0) = k°(q). The function g* (q reaches its 
minimum at the point e*^ q 0 ^, belonging for some n = n(q) to one of the intervals I n = I+U I~ (see the definition (l44l) l. 
Assume for instance that e*^ q = [^n+i> ^n\ ■ In this interval, the two dominant harmonics among the primary 

resonances are 5d = so(n) = Si and S 2 = so(n + 1). By the inequality (1481) . the second dominant harmonic among 
all resonances is S 2 = S 2 , at least for e £ 1+ close enough to e' n+ i, and we have to check that this is also true on 
the whole interval I+. Otherwise, assume that S 2 = s(q, 0) (a secondary resonance) for some values e £ Z+ (far 
from e' n+ i). Then, there would be an intersection between the graphs of g n+1 and g*^ q 0 ) in the interval I ++1 an dj 
view of the uniqueness given by Lemma ITT1 we would have g*^ q 0 ^(£ n ) < ~g n + i(£n) = Ji- A symmetric discussion can 
be done for the case £*^ q 0 j £ X~. 

By the above considerations, we have to check that, for any essential primitive q, and denoting n = n(q) as above, 
we have the lower bound 

s: ( 9 ,„)(e»)> Ja. (50) 

1/2 

Since the minimal value of the function g*^ q is ( 7 *) , it is enough to consider the essential primitives such that 

l /2 

( 7 *) < J 2 (by the lower bound (1301) . there is a finite number of such primitives). We have carried out a numerical 

verification of (15U1) for all metallic ratios with 1 < a < 10 4 . 
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Finally, we have to justify the the statement concerning the number of zeros 9* of the splitting function A4(9), for 
any e except for a small neighborhood of the transition values e. Notice that, since the second dominant harmonic 
changes from So(n — 1) to so(n + 1) as e goes from T~ to X+, the transition values are e = e n . As mentioned above, 
the dominant harmonics S\ = S\ = so(ro) and S 2 = S 2 = so(^A 1) are two consecutive resonant convergents, and 
hence we get k = 1 in (l62l) . As explained in Section H~2l this implies directly that the number of zeros of A4(9) is 
exactly 4k = 4. □ 

Remark. The exact equality in (15U1) holds true for some primitives q, as one can see in Figure []] for the concrete case 
of the bronze ratio. 


To end this section, we point out that the difficulty in this kind of results comes from the fact that we are considering 
infinite families of frequency ratios. But for a particular frequency ratio, it is always possible to obtain rigorous results, 
after a finite analysis of primary and secondary resonances. As an example, for Q = [ 1, 2, 2] it is not hard to see that 
the functions h\{e) and / 12 (e) have 4 and 12 corners in each period, respectively, as well as to give explicit formulas 
for their minimum and maximum values (this is illustraded in Figure [3]). 


4 Justification of the asymptotic estimates 

4.1 Approximation of the splitting potential by its dominant harmonics 

The last part of this paper is devoted to the proof of Theorem [I] which gives exponentially small asymptotic estimates 
of the maximal splitting distance and the transversality of the splitting. We start with describing our approach in a 
few words. 

Notice that Theorem |T] is stated in terms of the Fourier coefficients of the splitting function A 4 = VC introduced 
in m■ We write, for the splitting potential and function, 

£{9) = Y C k cos((fc, 9) — r fc ), M(9) = — Y Ah sin((/c, 9) — r fc ), (51) 

te 2 \{o} fce^\{o} 

with scalar positive coefficients C k , and vector coefficients 


Mk = kCk- (52) 

Although the Melnikov approximation @ is in principle valid for real 9 , it is standard to see that it can be extended 
to a complex strip of suitable width (see for instance |DGS04j '). from which one gets upper bounds for \C k — [iL k \ and 
| T k — cr k I, which imply the asymptotic estimates given below in Lemma [Till ensuring that the most dominant harmonics 
of the Melnikov potential L(9) are also dominant for the splitting potential C{9). The asymptotic estimates for the 
maximal splitting distance and the transversality, given in Theorem[I] are determined from a few (one or two) dominant 
harmonics of the potential. Thus, we consider approximations on C{9) given by such dominant harmonics, together 
with estimates of the sum of all other harmonics, which show that they are dominated by the most dominant ones. 

For the proof of part (a) of the theorem, that provides an asymptotic estimate for the maximal splitting distance, 
it will be enough to consider the approximation given by the first dominant harmonic. Thus, we write 

C(9) = C^\9)+^ 2 \9), £W(6)~/: Si cos((S u 9)-ts 1 ), (53) 

and we give below, in Lemma fl3l an estimate of the sum of all harmonics in the remainder J-^ 2 \9). This ensures that 
the maximal splitting distance can be approximated by the size of the coefficient of the dominant harmonic Si = Si (e) 
(see the proof of Theorem QJa) below). 

On the other hand, for the proof of parts (b) and (c) of TheoremQ] which concern the transversality of the splitting, 
we need to detect simple zeros of A 4(9). This is not possible with the approximation (1531) given by only one harmonic, 
and we need to consider at least two harmonics. Recalling that, by (l52l) . each (vector) harmonic A i k of M.{9) lies 
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in the direction of the integer vector k, we consider the two most dominant essential harmonics, given by linearly 
independent quasi-resonances S\ = Si(e) and S 2 = S 2 (e) (recall the definition of essential quasi-resonances at the 
beginning of Section [2.2l) . In fact, some non-essential harmonics c5i,c=2,...,m, can be (eventually) more dominant 
than S 2 - In order to show, in Section 14.21 that such non-essential harmonics have no effect on the transversality, we 
consider them separately, with specific upper bounds. We define the index of non-essentiality m = m(e ) > 1 as the 
integer satisfying 

ffSi( e ) < • • • < 9mSi( e ) - 9 s 2 ( £ ) < 9(m+l)Si- 

Recall from (H^l) that g* Si = cg* Si . It is clear that to = 1 if and only if the two most dominant harmonics are the 
non-essential ones S\ and S 2 - For instance, we see from Figure [ jf a) | that, for the case = [1,3], we have m = 2 for 
e belonging to some intervals, and m = 1 for the remaining values of e. 

Now, we write 


C(9) = C^(6) + T®{6) + T {3) (9), (54) 

C (2 \d) := £ Sl cos((Si,0) - t S i ) + Cs 2 cos((S 2 ,9) - rs 2 ), 

m 

T®(6) := E £ cSl cos(c(5i,0) - t cSi ), 

c—2 

and 7-^(9) containing all harmonics not in C^ 2 \9) or 7-^(9) (of course, we consider 7-^ = 0 if to = 1). Then, 
suitable estimates of the harmonics in T^ 2 \9) and 7-^(9), allow us to establish the existence of simple zeros 9* 
of 74(9), together with an asymptotic estimate for the minimal eigenvalue of D 74(9*), which can be taken as a 
measure for the transversality of the splitting (see Section m- However, we have to exclude some intervals where 
the second and third essential dominant harmonics are of the same magnitude and the approximation (1541) does not 
ensure transversality. Such intervals are small neighborhoods of the transition values e, where a change in the second 
dominant harmonic takes place. Such transition values can be defined as the values where 

h 2 (e) = / 13 (e). (55) 


We will use the following lemma, analogous to the one established in 1DG03 , DG04j . providing an asymptotic 
estimate for the dominant harmonics Cs^ and Cs 2 (and an upper bound for the difference of their phases ts, with 
respect to the original ones as*), as well as an estimate for the sum of all the harmonics in the remainders 
i = 2,2, 3, appearing in (l53l) and (l54l) . To unify the notation, we write ^ (•■■), defining the sets of indices 

k£Zi 

Z 2 = Z\{Q,S l ], Zs = {2S 1 ,---,mS 1 }, Z 3 = Z \ ({0, S u S 2 } U %). 

The estimate for each sum is given, due to the e^qronential smallness of the harmonics, in terms of the dominant 

harmonic in each set Zj, that we denote as S) = Si(e), i = 2,2,3. Notice that, for i = 2, the dominant harmonic 

is clearly Sj = 2 Si (non-essential) and, for i = 2, 3, the dominant harmonic can be either Si = Si (essential) or 

Si = 2Si-i (non-essential). With this in mind, we introduce the functions 


%(e) == 2/ii(e) =0|J £ )- 

hi(e) :=min (h z (e), 2h i _ 1 (e)) = 3 |.( £ ); * = 2,3, 


(56) 


We stress that, in the three cases, the function hi(e) is given by the minimum of the values g^.(e), with k belonging 
to the corresponding set of indices: hi(e) = mingle), i = 2,2,3. Comparing with the functions hi(e) defined 

k£Zi 

in (l46l) . we see that non-essential harmonics are also taken into account in the definition of hi(e). Notice also that the 
equality (1551) characterizing the transition values can be rewritten as h 2 (e) = /13(e). 


Recall that the coefficients Ct- introduced in m, are all positive. In fact, we are not directly interested in 
the splitting potential C(9), but rather in some of its derivatives (such as 74(9), T>74(9)). The constant Co in the 
exponentials is the one defined in (1391) . On the other hand, recall that the meaning of the notations and ' W has 
been introduced at the end of Section 11.31 
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Lemma 13 For e small enough and p = £ r with r > 3, one has: 


(a) C Si ~ p L Si ~ exp - 


C 0 hi(£) 

=- 1/4 


ks, - o- Si | ^ — 


= 1 , 2 ; 


keZi 


(b) > J £fc ~ -^4 % ~ 7174 ex P 




1/4 


Co^(e)] 
e 1 / 4 J’ 


* = 2,2,3. 


Sketch of the proof. We only give the main ideas of the proof, since it is similar to analogous results in I DG041 
Lemmas 4 and 5] and [DG031 Lemma 3]. At first order in p, the coefficients of the splitting potential can be approx¬ 
imated, neglecting the error term in the Melnikov approximation ©, by the coefficients of the Melnikov potential, 
given in (1571) : Ck ~ pLk = pctk e~^ k . As mentioned in Section 13.11 the main behavior of the coefficients Lk (e) 
is given by the exponents /3fc(e), which have been written in (1391) in terms of the functions <?&(£). In particular, the 
coefficients Ls t associated to the two essential dominant harmonics k = Si(e), i = 1,2, can be expressed in terms of 
the functions hi (e) introduced in (1461) . In this way, we obtain an estimate for the factor e _/3s <, which provides the 
exponential factor in (a). 

We also consider the factor a*,, with k = Si(e). Recalling from (H71) that |Sj| ~ £ -1 / 4 , we get from (1571) that 
a Si ~ £ -1 ' 4 , which provides the polynomial factor in part (a). 


The estimate obtained is valid for the dominant coefficient of the Melnikov potential L(9). To complete the proof 
of part (a), one has to show that an analogous estimate is also valid for the splitting potential C{9 ), i.e. when the 
error term in the Poincare-Melnikov approximation © is not neglected. This requires to obtain an upper bounds 
(provided in [DGS041 Th. 10]) for the corresponding coefficient of the error term in © and show that, in our singular 
case p = e r , it is also exponentially small and dominated by the main term in the approximation. This can be worked 
out straightforwardly as in |DG041 Lemma 5] (where the case of the golden number was considered), so we omit the 
details here. 


The proof of part (b) is carried out in similar terms. For the dominant harmonic k = Si inside each set Zj , 


i = 2,2, 3, we get 


£ 1 / 4 as in (S3- and an exponentially small estimate for Cz. with the function /q(e) defined 


in (156|l . Such estimates are also valid if one considers the whole sum in (b), since for any given e the terms of this sum 
can be bounded by a geometric series and, hence, it can be estimated by its dominant term (see |DG041 Lemma 4] for 
more details). □ 


In regard to the proof of Theorem ©a, c), we need to measure the size of each perturbation F <1 ' 1 (9) in (1531l54l) with 
respect to the coefficients of the approximations (9). Since by Lemma [T51 the size of F^ l \0) is given by the size of 
its dominant harmonic, we introduce the following small parameters: 


£Si 

Vi,j ■= 7-exp 

C Sj 


Co(hi(E)- 

£ l/i 



(/,j) = (2,1), (2,1), (3,1), (3,2), 


(57) 


as a measure of the perturbations in (l55H54]l . relatively to the size of the essential dominant coefficients Cs :j 
(we consider r/g 1 = 0 if the index of non-essentiality is m = 1). Although we define the parameters rjij in terms of the 
coefficients of C(9), we can also define them from the coefficients of its derivatives, such as the splitting function A 4(9), 


in view of (l52l) and the fact that the respective factors have the same magnitude: | S 3 \ 



The parameters r)ij always exponentially small in e, provided we exclude some small neighborhoods where Csj 
and Cg can have the same magnitude. For instance, we have 773 2 exponentially small if £ is not very close to the 
transition values (|55|) , at which the second and third dominant harmonics have the same magnitude. Analogously, the 
parameter 772,1 is exponentially small excluding neighborhoods where the first and second dominant harmonics have 
the same magnitude. 

Proof of Theorem\J\(a) . Applying Lemma IT51 we see that the coefficient of the dominant harmonic of the splitting 

function A 4(6) is greater than the sum of all other harmonics. More precisely, we have for £ —> 0 the estimate 

max |A4((9)| = |A4 Sl | (l + 0( V2 ,i)) ~ |AT Sl | ~ |Si|£ Sl , (58) 
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which implies the result, using the asymptotic estimate (S3 for |Si|, and the asymptotic estimate for I-Msj, in terms 
of hi(e), deduced from Lemma [Ha). 

We point out that the previous argument does not apply directly when e is close to a value where hi and h 2 
coincide, i.e. the first and second dominant harmonics have the same magnitude (for instance, for a metallic ratio Q 
this occurs near the values e' n , see fiH) and Figure [ jfb)] ). Eventually, more than two harmonics (but a finite number, 
according to the arguments given in Lemma UTl) might also have the same magnitude and become dominant. In such 
cases, the parameter 772,1 is not exponentially small, but we can replace the main term in (1581) by a finite number 
of terms, plus an exponentially small perturbation, and by the properties of Fourier expansions the maximum value 
of |A4(0)| can be compared to any of its dominant harmonics. □ 


4.2 Nondegenerate critical points and transversality 

This section is devoted to the study of the transversality of the homoclinic orbits for values of the perturbation 

parameter e, not very close to the transition values e defined in (l55l) . For such values of e we show that, under suitable 

conditions, the splitting potential C{9) has 4k nondegenerate critical points for some integer k > 1, i.e. the splitting 
function A i(9) has 4k simple zeros, which give rise to 4k transverse homoclinic orbits. Such critical points are easily 
detected in the approximation £ ( 2 '(0) introduced in (1541) . given by the 2 essential dominant harmonics, and using the 
estimates for J r ( 2 - ) (0) and J-( 3 \9) given in Lemma [T3l we can prove the persistence of the critical points in the whole 
function C(9). 

In fact, we make the computations easier by performing a linear change on T 2 , taking £^(9) to a very simple 
form. As in IDG031 IDG04j , we introduce the variables 

Vh = (Si,9) - r Sl , ip 2 = (S 2 ,9) - t S2 . (59) 

This change of variables is valid for e in the interval between two consecutive transition values, in which we have two 

concrete essential dominant harmonics S'i(e) and S 2 {e), which remain constant in this interval. In the new variables, 
the functions £, C^ 2 \ J 7 ®, J in (l5ll) become 

/C(t/>) = /C ( %0 + S (2 ) (V>r) + £ (3 ) (t/>) (60) 

where, in particular, we have 

/CW^) = £ S1 cos ipi + Cs 2 cos ip 2 . (61) 

It is clear that /C( 2 ) has the 4 critical points, all nondegenerate: ip *’° := (0,0), ( 7 r, 0), (0, 7 r), ( 7 r, 7 t) (one maximum, 
two saddles and one minimum, respectively). Regarding K,{ip) as a perturbation of /C ( 2 ^(t/>), we are going to show that 
it also has 4 critical points ip*, all nondegenerate, which are close to the critical points ip*’° of K,^ 2 \ip). We point out 
that, in general, the change (15^1) is not one-to-one on T 2 , but rather “K-to-one”, where 

k = k(s) := |det(Si, 5 2 )| ■ (62) 

Hence, the number of critical points of C{9) is 4k. It is not hard to show that n(e) is 4 In A-periodic in In e. Moreover, 
it is “piecewise-constant” with (eventual) jump discontinuities when changes in the dominant harmonics take place. 

Remark. For a metallic ratio f2 = [a], we know from Numerical Result [3(b) that k = 1 (a result checked numerically 
for 1 < a < 10 4 ), and hence there are exactly 4 transverse homoclinic orbits, for any e small enough (excluding a 
neighborhood of the transition values e). Although for other frequency ratios it is possible, in principle, to have k > 2, 
we have obtained k = 1 for all the cases we have explored. 


To establish the persistence of the critical points, we are going to use the following lemma, whose proof is a simple 
application of the 2-dimensional fixed point theorem and is omitted here. 
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Lemma 14 If /i ,/2 : T 2 


are differentiable and satisfy 

2 


fi + 


dfi 


dipi 


dfi 


dip 2 


< 1 , 


* = 1 , 2 , 


£/ien the system of equations 


sin ipi = fi(ip), sin ip 2 = / 2 WO 


(63) 


/ias exactly f solutions ip*, which are simple. Furthermore, if fi(ip), f 2 (ip) — 0 ( 77 ) f or an y V’ 6 T 2 , with 77 sufficiently 
small, then the solutions of the system satisfy ip* = ip*’° + 0{rf), with ip*’° = (0,0), ( 0 , 7 r), ( 7 r, 0 ), 


In order to apply this lemma, we introduce the following perturbation parameter, using the parameters rjij defined 
in (EZD, 


V ■= max(77 2 l , 773 , 1 , 773 , 2 ). 


(64) 


Lemma 15 The function K,(ip) has exactly f critical points, all nondegenerate: 

ip* = ip*’° + 0(ri), with ip*’° = ( 0,0), ( 0 , 7 r), ( 7 r, 0 ), ( 7 T, 7 r). 
At each critical point, we have 


(65) 


det DK.(ip*) = S15* 2 C Sl Cs 2 (1 + 0( 77 )) 


where 5* = cosip*’ 0 = ± 1 , i = 1 , 2 . 


Proof. We see from ((60][6T|) that the system of equations S7K.(ip) = 0 can be written as in (|63|) . with the functions 


A (VO 


1 ( d g® dg^ \ 

y dipi + dip! J ’ 


AWO 


1 < 9 £ (3) 

Cs 2 dip 2 


(notice that does not depend on ip 2 ). By Lemma fTTTl we have fi(ip), f 2 {ip) = 0 ( 77 ), with 77 as given in (l64ll . Hence, 
applying Lemma ITT1 we deduce the result for the critical points of K.(ip). 


We also provide, for each critical point, an asymptotic estimate for the determinant of D 2 K.{ip*). It is clear, for 
the perturbed critical points, that the signs 5* = ±1 become perturbed as follows: cos ip* = 6* + Off 2 ). Writing 

o / fcn k\ 2 \ _ 

D IC(ip) = , we have from (1601) the approximations 

K22 / 


4-12 


fell = 
hi2 = 
^22 = 


d 2 IC 

W 2 


1 

d 2 K. 


-£ Sl ( COSIp! + £>( 77 ^, 773 , 1 )) , 


= £-Si ■ 0(773,1) = Cs 2 ■ 0(773,2), 


dipi dip 2 
d 2 K. 

—3 = -Cs 2 (cos 1 p 2 + 0 ( 773 , 2 )) 
dipf 


and we deduce the expression for the determinant of D 2 /C(t/>*). 


( 66 ) 

(67) 

( 68 ) 

□ 


Now we complete the proof of part (b) of our main theorem by applying the inverse of the linear change (|59l) to 
the critical points ip* of IC(ip), in order to get the critical points 0* of the splitting potential C(9), i.e. the zeros of the 
Melnikov function AA(9). 


Proof of Theorem \J\(b). Since the linear change (15U1) is “fc-to-one”, with k as in (1(721) . the 4 critical points ip* 
of IC(ip) give rise to 4 n critical points 9* of C{9). It is clear that such critical points are also nondegenerate, and hence 
they are simple zeros of the splitting function M.(9). □ 


Remarks. 
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1. Recall that the vectors Si = Si(e) remain constant between consecutive transition values e (see the proof of 
Proposition [12]). On the other hand, by © the points ip* are 0{rj)- close to the points tp*’°, where rj is 
exponentially small. Hence, the points 9* = 9*(e) remain “nearly constant” along each interval of e between 
consecutive transition values e, and can “change” when e goes across a value e. 

2. As a particular interesting case, we may consider the phases tJk = 0 in the perturbation ©. In this case, our 
Hamiltonian system given by <(4] [9]) is reversible with respect to the involution 


K (~x,y,-(p,I) (69) 

(indeed, its associated Hamiltonian field satisfies the identity Xh o 1Z = —TZXjj)- We point out that reversible 
perturbations have also been considered in some related papers |Gal94l IGGM99b l I RW98 j. Under the reversibil¬ 
ity (l69l) . the whiskers are related by the involution: W s = 1ZW U . Hence, their parameterizations in (1121) can 
be chosen in such a way that J s (6) = J u (—9 ), provided the transverse section x = tt is considered in their 
definition. This implies that the splitting function is an odd function: Xi(—9) = —A4(9) (and the split¬ 
ting potential C{6) is even) and, using its periodicity, one sees that M{9) has, at least, the following 4 zeros: 
9* = (0,0), (7r, 0), (0,7r), (7 t, 7t) (notice that they do not depend on e). Although such zeros could be non¬ 
simple in principle, the result of TheoremQJb) says that they are simple for any e except for a small neighborhood 
of the transition values e, at which some bifurcations of the zeros could take place. 


It remains to provide, for each zero 8* of the splitting function A4(0), an asymptotic estimate for the minimal 
eigenvalue of the matrix D 2 C(9*) = DA i(8*), as a measure for the transversality of the splitting. 


Proof of Theorem\l\(c) . Denoting D = detD 2 C{9*) and T = trD 2 £(#*), we can present the (modulus of the) 
minimal eigenvalue of D 2 C{9*) in the form 


TO = 


2|D| 


\T\ + \/ T 2 - ID 


\D[ 

\T\ ’ 


where we have taken into account that 0 < \JT 2 — 4 D < \T\. Thus, we need to find estimates for |D| and |T|, at 
the critical points 9* of C(9) (or zeros of A4(0)). 


By the linear change (|59l) . we have D 2 C{9*) = A T D 2 JC(ip*) A , where A is the matrix having the vectors Si and S 2 
as rows. In (|62D . we have defined k = |det A\. Since k = k(e) is picewise-constant and periodic in lne, it is bounded 
from below and from above: n ~ 1. Applying Lemma fl5l we get the asymptotic estimate 

\D\ = Cs 1 Cs 2 (1 + 0(r])) — C Sl C S2 . 


On the other hand, in order to estimate T we write D 2 /C as in (l66l l68l) . and obtain 

D 2 C{8) = k n Si ■ Sj + ki 2 (Si ■ Sj + S 2 -Sj) + k 22 S 2 ■ Sj, 

which implies that T = kn (S'lll + 2ki 2 ( Si,S 2 ) + k 2 2 IS 2 I 2 , where |-| 2 denotes the usual Euclidean norm (which is 
equivalent to the norm |-| = 1-^ mainly used in this paper). Now we use, at the critical points ip*, the estimates for the 
matrix D 2 IC(ip*) given in (l66ll68ll . We obtain |/cn| ~ Cs 1 as the main entry, and |A 12 I ~ Cs 1 ■ 0(rj 3 ^), \k 22 \ ~ Cs 2 - 
Applying also the estimate (l47l) . we obtain 


\T\~-=C Sl , 

v £ 


and hence to* 


\D[ 

\T\ 


y/~E Cs 2 ■ 


Applying the estimate for Cs 2 given in Lemma fT51 we obtain the desired estimate for the minimal eigenvalue. 


□ 
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